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ABSTRACT 

The purpose of this thesis is to develop the equations of condi- 
tion necessary for determining the coefficients for Runge-Kutta methods used in 
the solution of ordinary differential equations. The equations of condition 
are developed for Runge-Kutta methods of order four through order nine. 

Once developed, these equations are used in a comparison of the local trun- 
cation errors for several sets of Runge-Kutta coefficients for methods of 
order three up through methods of order eight. 
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CHAPTER I 


INTRODUCTION 

1.1 Description of the Problem 

i 

“The purpose of this thesis was to develop the equations of condi- 
tion necessary to determine the coefficients for Runge-Kutta methods used in 
the solution of ordinary differential equations. The equations of condition 
were developed for Runge-Kutta methods of order four through order nine. 
Once, developed , the equations were used in a comparison of local truncation 
errors for several sets of Runge-Kutta coefficients for methods of order 
three through order eight. The equations of condition were generated by 
the computer using the algebraic manipulation language SYMBAL 3 and the 
CDC 6600/6400 computer system at the University of Texas at Austin. 

1.2 Motivation 

Numerical solutions to ordinary differential equations play a 
large role in the study of science and engineering. The integration of 
spacecraft trajectories in orbital mechanics, for example, requires ac- 
curate and efficient numerical integration techniques. The Runge-Kutta 
type methods for solving systems of ordinary differential equations are 
widely used because of their simplicity of application. Any studies of 
existing Runge-Kutta methods , or attempts to produce new or better Runge- 
Kutta methods will result in a better understanding of the methods, and 
will lead to the development of more efficient techniques for solving the 
ordinary differential equations of science and engineering. 
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CHAPTER II 


GENERATION OF THE EQUATIONS OF CONDITION FOR RUNGE-KUTTA 
COEFFICIENTS WITH TAYLOR SERIES EXPANSIONS 


2.1 Description of the Method 

Consider the ordinary differential equation 

= y'(x) = f(x, y) (1) 

with initial conditions y(x Q ) = y , where y and f may be vectors, y is the 
dependent variable, and x is the independent variable. A solution of the 
form y = y(x) which satisfies the initial conditions and Equation (1) is 
desired. The existence and uniqueness of a solution is assumed since there 
exist theorems 8 which state that if f(x, y) is sufficiently well behaved 
near a point (x, y), then Equation (1) has a solution that passes through 
the point and is unique in a neighborhood of the point. 

A solution of Equation (1) can be found using a Taylor series ex- 
pansion of y about y = y^ in the form 


y(x + Yh) = y + hyy ’ + 


2 ., 2 


h z y 


+ nfr 

' ~ x or 


"l" • • • • 


( 2 ) 


o J o ■ " ,J o 21 J o 3! ^o 

for any y for which the series converges. This technique requires, however, 
that the derivatives of y(x) be known up to some desired order. Since 
f(x, y) is a function of both x and y, these derivatives may become quite 
complicated. It is convenient to introduce the following notation: 

y’ = = f , 

y" = f^- 1 -* 
v" * = f*- 2 -' 


2 
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where ^ _ + 5 

where f and f represent partial derivatives of f(x, y) with respect to 
x y 

x and y respectively. Successive derivatives of y(x) can be generated using 

f [j + l] = f Cj] + 

x 

If the function f(x, y) is complicated, this procedure can become prohib- 
itively difficult. 

Runge* 7 was the first to point out that it was possible to avoid 

the successive differentiation required in the Taylor series solution while 

still preserving the accuracy. Runge bypassed the derivatives in the Taylor 

series solution by using evaluations of the function f(x, y) within the 

interval (x , y ) to (x + h, y(x + h)). Runge 's ideas were applied to 
o o o o 

first order differential equations in a more accurate form by Heun and 
Kutta*®, and extended to second order differential equations by Nystrom . 
Zurmuhl 20 continued the extension to nth order. 

Considering only the first order system of differential equations, 

the problem formulation now becomes 

y' ■ f(x> y) ■ 

y'V = v o 

with f Q = fU 0 » y o )j and 

k-1 

f k = f(x o + °k h ’ y o + 6 kX f X } » k = I* 2, 3, .... n, 

where n + 1 equals the number of function evaluations required in the inter- 
val. The solution is then given by 

Y(x) = y 0 + h l c^f^ + 0(h m+1 ) 


( 3 ) 



where m is the order to which the Runge-Kutta formula agrees with the Taylor 
series, and h is the integration step-size which implies a value y = 1 in 
Equation (2). 

In order to obtain Y(x) from the Runge-Kutta scheme, it is neces- 
sary to determine the coefficients a, 3, and c in such a way that the 
Runge-Kutta solution is equivalent to the Taylor series solution to some 
order m. In order to accomplish this, the common approach in the past has 
been to equate the expression for y(x Q + h) obtained from the Taylor series 
in Equation (2) and the expression for Y(x) obtained from the Runge-Kutta 
formulation given in Equation (3). From this equivalence the coefficients 
of corresponding powers of h are compared yielding a system of nonlinear 
algebraic equations, which are referred to as equations of condition, from 
which the unknown coefficients may be determined for the classical Runge- 
Kutta formula of order m. 

For the solution of this set of equations of condition to yield 
coefficients which produce a method of order m, there is a minimum number 
of function evaluations required for that order. For example, since a 
seventh order method requires at least nine function evaluations , a method 
which has eight function evaluations can at best be a sixth order method 
although only seven function evaluations are actually required for a sixth 
order method. Since many of the existing methods use more than the minimum 
number of function evaluations for a given order, the number of function 
evaluations used for each order in this study was chosen equal to that of 
the method using the most function evaluations for that order. Table 1 
shows the minimum number of function evaluations , the number of function 
evaluations used in this study, and the number of new equations of condi- 
tion expected for methods of order three through order nine. 
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Table 1. 

Number of Function Evaluations and Equations of Condition 


Order 

Min. No. of Func. 
Eval . Reqd . 

No. of Func. 
Eval. Used 

No . New Eqs . 
Expected 

3 

3 

4 

2 

4 

4 

5 

4 

5 

6 

' 6 

9 

6 

7 

8 

20 

7 

9 

11 

48 

8 

11 

13 

115 

9 

— 

17 

286 


It will be advantageous to give an example to illustrate this 
procedure, which, although basically simple and straightforward, is ex- 
tremely tedious . Consider the development of the equations of condition 
for a Runge-Kutta formula of third order with three function evaluations. 

The Taylor series solution becomes : 


y (x 0 

+ h) = y o + hy; 

+ ~ y" 
2 y o 

+ b y o' + 0(h4) * 

where 

= f o ’ 



y 1 f 

= (f + f f) 



J o 

x y o 



v T 1 ' 

= (f + 2ff 

+ f 2 f 

+ f f + ff 2 ) . 

J o 

xx xy 

yy 

x y y O 


By expanding the function f^ in a Taylor series of two variables, the 
Runge-Kutta solution yields: 

Y(x) = y Q + h ( c Q f Q + c 1 f 1 + c 2 f 2 ) + 0(h 4 ) , 

where f Q = f Q , 
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f O + “2 hf X + h[6 20 f * S 21 <f + “l hf X + h6 10 f y f) ] f y 


* a 2 h2[B 20 f + e 21 (f o * “l hf x )]f xy * 0<h3) - 


f i = f< *o + a i h ’ y 0 * he io f o> 

' f o * “l hf x + hB 10 f y f * ^l h2f xx + ^lO^^yy + 0<t ‘ 3) - 
f 2 * f(x c, * “2 h - ^ + h<S 20 f c + B 21 f l )) 

= f o + “2 hf x + h [ e 20 f + 6 21 (f + “l hf x + h6 10 f y f, ] f y 

* V 2 h2f xx * *> 2 l>20 f o + ®21 (f o * "l hf x'] 2f y, 

* a 2 h2[B 20 f + e 21 (f o * “l hf x )]f xy * ° <h3 )' 

Setting y(x Q + h) equal to Y(x) gives: 

y o + hf o + h2 [|< f x + f y f)] + h3t f <f xx + 2f xy f + f2f yy * f x f y + ff y 2>:i * 0<h ‘' ) 

= y o + h{c o f o + c l [f o + h S f x + 6 10 f y f) * ^ 2( “i f xx + *i 0 f 2 V 

+ h< “l 6 10 ff xy ): > + c 2 [ f o + h(o 2 f x + S 20 ff y + 6 21 f y (f + “l hf x 

+ S 10 f y ft)) + 7 h2<a 2 f xx + f2f yy (6 20 + ®21 ,2) + h2(o, 2 6 20 ff xy 
+ OjB^ff )]} + 0(h4) ‘ 

Comparing similar powers. of h: 

v = v : 


y c = 

f = 

y o ; 
c f + 

c, f + c 

o 

O O 

1 o 

¥ = 

c, a, f 

+ c^a^f 

2 x 

1 1 X 

2 2 x 

■f f = 

! y 

C l 3 10 f y 

f + c 2 ($, 


Kx = 7 (a i=l + a 2 c 2 >f xx • 

^ 2f yy = 7 (c l 6 l5 + c 2 (6 20 + 6 21 )2)f2f yy . 

K f y = c 2 B 2l“l f x f y • 

H* = C 2 6 21 6 10 f y f • 

Ky f = [c 2“2 (6 20 + *21 J + c l“l 6 10 ]f xy f • 
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An examination of these equations reveals that 


a l = 6 10 • and 
°2 " B 20 + e 21 * 

In general, this requirement can be stated 

i-1 

a. = 7 0 . . . 

1 j=0 ^ 

The equations of condition for the third order Runge-Kutta formula then 
become: ; 

. ' 1 = c o + C 1 + c 2 ' 


2 = °l“l + °2 a 2 

t ~ c i a i + ° 2 “2 


6 " C 2 a l 6 21 


By extending this procedure to higher orders it is observed that 
as the order and number of function evaluations increase, the expressions 
^ Taylor series and f^ in the Runge-Kutta formula become more 

and more complex and grow enormously in length. 

2.2 Computational Procedure and Limitations 

In order to generate the equations of condition for Runge-Kutta 
methods, the computer and the algebraic manipulation language, SYMBAL, were 
employed. From the example just given it is observed that the generation 
of the equations of condition for even the fourth and fifth order methods 
is an extremely tedious task. Therefore, the computational power of the 
computer was used to develop these equations, beginning with the fourth 
order method, in hopes of eventually developing the equations of condition 
for high order Runge-Kutta methods . 
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The SYMBAL (SYMBolic ALgebra) language used was developed by 
M. E. Engeli 4 and is implemented on the CDC 6600 at The University of Texas 
at Austin and at the Swiss Federal Institute of Technology in Zurich. The 
language is a generalization of ALGOL 60 with the manipulation of unre- 
stricted algebraic expressions of which numbers represent only a very 
special case. Although the language SYMBAL may not be widely available, 
it. is similar to several other algebraic manipulation languages in use 
which may be more familiar to the reader. 

The original program used in this study is shown in Appendix I 
and appears in the SYMBAL manual 4 where, as an example, the fourth order 
equations of condition are developed. This program proceeds exactly as 
in the third order example of the previous section. As a result, the 
Taylor series expansions are in two variables and the requirement that 
a = SB was not included. The first revision of this program. Revision 1 
shown in Appendix I, incorporates an. autonomous differential equation to 
avoid the Taylor series expansions in two variables, i.e., the right-hand 
side of the differential equations contain only the dependent variable. 

The fact that a = Eg was also incorporated into this program. Table 2 in 
the next section gives the results of this change. 

With the increased efficiency gained by the first revision for 
the fourth order Runge-Kutta equations of condition, efforts were then 
directed toward generating the fifth order equations of condition. The 
fifth order method requires six function evaluations and all expansions 
were carried out to. include fifth order terms , thus increasing the complex- 
ity of the problem considerably over that of the fourth order. With the 
increased complexity of computation, problems of inefficiency in the pro- 
gram became apparent. Due to the size of the expressions and to the large 
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number of intermediate calculations, storage problems were encountered. 

When it became evident that this would be a serious problem, the program 
was divided into two parts: (1) y(x Q + h) and Y(x) were calculated as 

described earlier, and (2) like powers of h in the two resulting expres- 
sions were compared to yield the equations of condition. 

The major area of inefficiency in the program was found to be 
the calculation of unnecessary terms in many instances which, after multi- 
plications or substitutions, led to expressions containing higher order 
terms in h than were needed. The program at this stage was set up to 
develop all needed expressions to fifth order, then multiplications and 
substitutions were carried out and the resulting expressions truncated to 
fifth order. A considerable effort was made to minimize the calculation 
of these extra terms so that only those needed to make the final result 
fifth order were retained. The details of this work will not be discussed, 
but a copy of this program. Revision 2, is shown in Appendix I. Attempts 
at further revisions of this basic program did not lead to substantial im- 
provement. The results of generating the equations of condition for the 
fifth order are presented in Table 3 of the next section. 

From the results shown in Table 3 it is obvious that this method 
would not produce the fifth order equations of condition. As a final at- 
tempt to determine the fifth order results with this method, efforts were 
made to use magnetic tapes to help decrease storage and to store expressions 
for future use. However, the SYMBAL compiler has very limited file manipu- 
lation capabilities. It was possible to store results on tape, but there 
are no provisions for reading information into a SYMBAL program from tape 
so that the usefulness of this effort was considerably reduced. 



JO 

Because this method was using the maximum allowable storage and 
because relatively long computing times were encountered at the fifth order 
level, it was evident that the series expansion approach was not practical, 
especially since the equations of condition were sought for high order 
formulas (eighth, ninth, tenth, etc.). Although this method failed to 
achieve the desired goal, it did provide much insight into the capabilities 
and limitations of the SYMBAL language. 

2.3 Summary of Results 

The results of the fourth order programs are given in Table 2. 

These results are for a complete set of equations with four function evalu- 
ations. Table 3 gives the results for the fifth order Runge-Kutta method 
with six function evaluations. It should be remembered that the fifth 
order program was divided into two sections, and Table 3 gives only the 
results for the first half of the program which was to calculate y(x Q + h) 
and Y(x). Although these expressions were finally obtained and stored on 
tape, the actual set of equations of condition was not obtained. The pro- 
gram, Revision 1, used in the fifth order study was the same as that used 
in the fourth order study with only the order and number of function eval- 
uations changed. 
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Table 2. 


Fourth Order Results 


Program 

CDC 6600 . 

No, of Eqs. 

Storage ReqcL 

Run Time 

Original 

77000 

8 

30.4 sec. 

19 

Revision 1 

77000. 

8 

10.8 sec. 

7 


Table 3. 

Fifth Order Results 


Program 

CDC 6600 

y(x Q + h), Y(x) 

Storage Reqd. 

Run Time 

Revision 1 

220000g (max.) 

500 + sec. 

not obtained 

Revision 2 

220000g (max.) 

97 sec. 

obtained 











CHAPTER III 


GENERATION OF THE EQUATIONS OF CONDITION FOR RUNGE-KUTTA 
COEFFICIENTS WITH THE METHOD OF E. BAYLIS SHANKS 


After the failure of the Taylor series approach discussed in 
Chapter II, an alternative method was sought. Since most of the problems 
in the Taylor series approach stemmed from the large expansions of f^^ 
and f^ it was natural to turn to the work of E. Baylis Shanks 19 . Shanks 
has developed a method which generates the desired equations of condition 
without carrying out the classical expansions and comparing coefficients 
of similar powers of h, as was done in Chapter II. 

3.1 Description of the Method 

Again the system under consideration is 

y* = f (*» y) » y(* 0 ) = y Q , 

where the solution is given by 


n 

Y(x) = y + h 7 c . f . 

i=0 1 


S 


with 


i-1 


f . = f<* o + . y o + h l e f ). 

3=0 J J 


The parameters to be determined are again a, g, and c. 

It can be shown 19 that a necessary and sufficient condition for 
the Taylor series for y(x Q + h) and the expansion of the solution Y(x) to 
agree through terms in h m is 


(f 


[k-1]. 


} o = k l c i (f i Lk “ 1J ) 0 ’ k = 1. 2 m-1, 

i=0 1 ° 


[k-1] 


(4) 


where f was defined in Chapter II, and f . ^ ^ denotes the kth derivative 

i-1 1 

of f L as expanded about (x q + a^, y + J g. .f.). Shanks states and proves 

j=0 ^ ] 


12 



13 


several theorems in order to arrive at expressions for the derivatives 

(f [k-i] > and (f.t k_1 ^) in Equation (4). He shows that the left-hand 
O 1 o 

side of Equation (4) may be represented by 


where the factor is a product of integers and partial derivatives. 

fkl 

For an expression for q Shanks first: defines a quantity 
*4 ^ m and he then develops inductively that the general term, Q. ^ 


fid s E tt ^ 
t t 


( 5 ) 


t 


IS 


aj!( nVj"B ..Q . Lk " 1J ) s 

1 r=l j=l 11 3 t r=l j=l 1] 3 t 


[ 1 ] 


A k, 


[k-1] 


( 6 ) 


where the induction is on k 3 k = i + 2 . + •••+- kd » an< ^ w bere the sub 

1 2 k 

scripts t, t are used to number distinct terms defined for each fixed k. 

As an example consider 1$ = 2 where either: 

il * 2, i 2 = 0, t = l, 

[ 2 ] 2 

q 4 - a i ’ 


or 


i 1 = 0, i 2 = 1, t = 2, 

q.M = 2 Ye.*^ • 

x 2 jSi 11 1 

It is then stated and proven 19 that fjj k ^ can be given as 

f[ k ^ = E .Q.^ k kp^ . 
i n t t 

From Equations (4), (5), and (7) the necessary and sufficient condition 


( 7 ) 


= (k + 1) f c i z t Q i [k] ^ k] . 


i=0 


can now be written: 
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or, by rearranging this expression, 

Z [l - (k + 1) l c Q = 0 . (ft) 

i=l t 

From Equation (8) the sufficient condition then becomes: 

n Tkl 

1 = (k + 1) l C.Q. LKJ , k = 1, 2,..., ift-1. (9) 

i=l \ 


The reader is referred to the original paper by Shanks 1 ^ for a 
complete treatment of the theorems and proofs leading to Equation (9). 
Equation (9) gives the sufficient conditions to be satisfied for a given 

Tkl 

order and these equations can be completely determined when the Q. J ex- 

X t 

pressions are known. 

For comparison with the Taylor series approach, the third order 
equations will be developed using Shanks' method. From Equation (4), 
with k = 1, it is seen that: 


but since f. = f 
i o 
o 



1 


n 



This equation always occurs and is thus assumed for simplicity from this 
point on since it is not produced by Equation (9). From Equations (6) 
and ( 9 ) , for k = 1 : 


Q W = a. 


1 = 2 J c.a. ; 
i=l 1 1 


and for k = 2 , 




1=37 c.a? , 
. L n 1 1 ’ 

i=l 


Q if 3 ■ 2 • 

1 = 3 j 2 °i 2 ji 6 n“i • 



The equations of condition are then: 


% + C 1 + C 2 • 


2 = C l°l + C 2°2 * 

■3 = C l a l + c 2 a 2 ’ 


6 C 2 B 2l a i 


which are identical to the equations obtained in the example in Chapter II. 
The series expansions have been eliminated and the computational complexity 
greatly reduced. /' 

3.2 Computational Procedure and Limitations 

A SYMBAL program was written to produce the Q expressions, from 
which the equations of condition follow immediately. This program is shown 
in Appendix I as Shanks' Method and includes the output for the fifth order 
equations of condition. It was found that Shanks' method was particularly 
well suited for computer implementation. 

Equations of condition were generated for the fourth, fifth, 
sixth, and seventh order Runge-Kutta methods without difficulty. The 
number of function evaluations used in each order is given in Table 1 in 
section 2.1. The results of this program are given in Table 4 in section 
3.4. From Table 4 the time required for successive orders is seen to grow 
rapidly. The time that would be required for an eighth order run was 
estimated at about eight minutes. However, it was not possible to run this 
program because of what appeared to be a limit on the size of an array or 
the number of variables used. The program could not be compiled, and was 
terminated during the setting up of the Q array which is a triangular 
array with three indices of the form Q [ 8 , 13, 297]. An "inventory over- 
flow" diagnostic was given and all efforts to compile the eighth order 
program failed. 
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3. 3 Generation of Reduced Systems of Equations of Condition for 
High Orders 

Up to this point the goal had been the generation of the complete 
set of equations of condition for each order of Runge-Kutta method. 

However, to solve the resulting equations for the a, 6, and c coefficients, 
the number of equations must be reduced for the higher order methods since 
the system is overdetermined, i.e., more equations than unknowns. The 
usual procedure in solving the equations of condition is to generate the 
complete set of equations and then make certain assumptions to reduce the 
number of equations . The assumptions that are usually made arise naturally 
from a quadrature approach 16 and they make several of the original equa- 
tions identical, thus reducing the number of equations in the system. 

The assumptions used in this study were: 

EBa = ia 2 , 

EBa 2 = -|a 3 , 

EBa 3 = -ia 4 . (10) 

These assumptions were incorporated into the calculation of the Q expres- 
sions which then led to the generation of a reduced system of equations 
directly. This program is shown in Appendix I as Shanks' Revision 1. 

The assumptions that reduce the number of equations of condition, 
and the number of these assumptions to be made to minimize the number of 
function evaluations required are discussed by Curtis 3 , who concludes that 
four assumptions of the type given by Equation (10) will produce an eighth 
order method with eleven function evaluations . Since the present study 
is concerned with generating the equations of condition, it was decided 
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to use only the three assumptions given in Equation (10) solely on the 
basis of ease of computation of the eighth and ninth order equations of 
condition. The eighth and ninth order equations of condition were then 
generated in reduced form without difficulty. The results are given in 
Table 5 in section 3.4. 

With this technique it is now feasible to go to the tenth and 
higher order methods by making additional assumptions, but this was not 
done since it was desired to use the equations generated up to this 
point to make comparisons of local truncation errors. This comparison is 
discussed in Chapter IV. It is stressed, however, that the equations of 
condition can be computer generated for higher order Runge-Kutta methods. 

3.4 Summary of Results 

The results for the fourth, fifth, sixth, seventh, and eighth 
order methods using the Shanks method are given in Table 4. These are 
complete sets of equations, but the number of equations given for each 
order in the table is the number of new equations due only to that order. 
Table 5 gives the results for the reduced eighth and ninth order sets of 
equations of condition. The number of equations given in Tables 4 and 5 
do not always agree with the number of equations expected for the higher 
orders since the SYMBAL program generates some identical equations for the 
higher orders. 



Table 4. 


Results from Shanks ' Method 



CDC 6600 


Storage Reqd. 


77000, 


77000, 


77000, 


150000, 


220000, 


Run Time 


1.4 sec. 


3.8 sec. 


12.0 sec. 


75.5 sec. 


Not Compiled 


No. of New Eqs. 




Table 5 . 


Results from Revised Shanks' Method 


CDC 6600 


Storage Reqd. 


77000, 


Run Time 


12.0 sec. 


120000, 


No\ of New Eqs . 














CHAPTER IV 


CALCULATION AND COMPARISON OF LOCAL TRUNCATION 
ERROR COEFFICIENTS 

4 . 1 Description of the Computational Procedure 

As a means of comparing sets of coefficients for Runge-Kutta methods 
of a given order, the coefficients of the local truncation error terms can 
be evaluated. From the development of the equations of condition, it is 
observed that these equations arise as multiplicative factors of the par- 
tial derivatives that appear in the expansions of the solutions y(x Q + h) 
and Y(x), as given in Chapter II. For these two solutions to agree through 
some order m, the equations of condition arising from comparing terms in h 
up through order m must be satisfied exactly. The leading term of the 
truncation error consists of the partial derivatives multiplied by numer- 
ical factors. These factors are simply the equations of condition that 

m+ j_ 

arise from the h terms in the expansions , and are explicit functions of 
the coefficients of the mth order Runge-Kutta method. More explicitly, 
these truncation error terms are of the form 

TE(i)[. . .] i h m+1 , 

where TE(i) represents the truncation error coefficient and [...] represents 
the partial derivatives which are multiplied by that particular TE(i) . The 
index i numbers the distinct terms in the h portion of the expansion. 

By evaluating these multiplication factors, or truncation error coefficients, 
and by comparing the results with those obtained from using other sets of 
Runge-Kutta coefficients for the same order method, an "optimum" set of 
coefficients for use with a given order Runge-Kutta method can be deter- 
mined by determining the set that has the lowest values for these trunca- 
tion error coefficients. 
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The equations of condition were generated in complete form for 
orders four through seven. These equations were used to produce the 
truncation error coefficients for the third through the sixth order methods. 
For the eighth and ninth order methods, it was necessary to generate re- 
duced sets of equations of condition, and these were then used to obtain 
the truncation error coefficients for the seventh and eighth order methods. 

The truncation error coefficients for several sets of Runge-Kutta 
coefficients are given for methods of order three through methods of order 
eight in Table 7 through Table 12 in the next section. 

4. 2 Comparison of Truncation Error Coefficients 

Several sets of Runge-Kutta coefficients for methods of order 
three through eight were used in this comparison. These coefficients are 
given in Appendix II for reference. Several points concerning notation 
should be emphasized for these sets of coefficients. The Fehlberg 5 ’ 6 
methods are given as "Fehlberg m(m + 1)" where m is the order of the method, 
with a solution of order m + 1 also being calculated to use with an auto- 
matic step-size control. The difference in the two solutions is used as 
an estimate of the truncation error made during that step. The new step- 
size is then based on this error and the desired accuracy. When this step- 
size control is used, the Fehlberg methods require extra function evalu- 
ations to produce the solution to both orders. The advantage of the auto- 
matic step-size, however, makes these methods desirable, especially if the 
system of differential equations is such that a rapidly varying step-size 
is required. Since these methods can be used to give a solution of order 
m or of order m + 1, they are shown in the comparisons for both orders. 

Table 6 shows the Fehlberg methods with the number of function evaluations 


o 
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required for each order solution with and without the automatic step-size 
control. 

The notation for the Butcher 2 and Shanks 19 methods are of the 
form ’’Shanks (m - n)” where m denotes the order of the method, and n de- 
notes the number of function evaluations required per step. It should be 
noted also that some of these methods by Butcher and Shanks are only ap- 
proximately of order m. They were developed such that the number of func- 
tion evaluations used is less than that normally needed for a true mth 
order method. For this reason the Shanks (6-6) is compared as both a 
fifth and a sixth order method, and the Shanks (7-7) is compared as both 
a sixth and a seventh order method. Methods which are approximations of 
this type are denoted with an asterisk in the tables. 

A large number of fifth order methods fall within a family of 
methods known as the Newton-Cotes Family. Only three of these forms are 
given although several others exist. The "UT" coefficients 1 are sets of 
coefficients developed at The University of Texas at Austin along the 
Fehlberg 4(5) format so they appear in both the fourth and fifth order 
comparisons. 

The truncation error coefficients are denoted by TE(i) and are 
given in Table 7 through Table 12 for the third through eighth order 
methods . The first row in each table gives the number of function eval- 
uations required (F.E.R.) per step for each set of coefficients. 
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Table 6 . 

Function Evaluations Required 
for the Fehlberg Methods 


Method 

Order 

of Solution 

j 

j No. of Function Evaluations Reqd. 

With Step-Size Control 

Without Step-Size Control 

Fehlberg 4(5) 

4 

6 

5 

Fehlberg 4(5) 

5 

6 

6 

Fehlberg 5(6) 

5 

8 

6 

Fehlberg 5(6) 

6 

8 

8 

Fehlberg 6(7) 

6 

10 

8 

Fehlberg 6(7) 

7 

10 

10 

Fehlberg 7(8) 

7 

13 

11 

Fehlberg 7(8) 

8 

13 

13 

Fehlberg 8(9) 

8 

17 

15 







TAdlE 7 

COMPARISON OF THIRD ORDER METROS 



FEHI.RERO 3(4) 

CLASSICAL 

RALSTON OPTIMUM 

HEUN 

NYSTROM 

E • E • R • 

4 

3 

7 

3 

3 

TE( 1 ) 

. '-r.233Q 

0,000000 

.02^833 

. 027777 

.027777 

r E ( 2 ) 

. 0.-2339 

. 0*3333 

0,000000 

.027777 

,027777 

T E < 3 ) 

0 „ o n n o o ft 

G.OOUOOO 

o ,006000 

.083333 

.083333 

f E ( 4 ) 

. ^263 1 S 

,250000 

,250000 

.250000 

,250000 



TABLF 8 


COMPaRISCN Of FOURTH ORDER VEthODS 


F.E 

• R • 

FEHI.PERr, 3(4, 
5 

FEHLeEHG U5)-l 
5 

FEHLRFRG 4 (5) -2 
5 

CLASSICAL 

4 

TE( 

i) 

.000249 

.000057 

.000020 

,000347 

f E ( 

2> 

,001498 

.000347 

,000120 

.002083 

te< 

3) 

.001388 

.000231 

,000080 

.002083 

TE ( 

4) 

.005774 

.002083 

,001782 

.008333 

TE ( 

5) 

. 000999 

.000231 

.000000 

.001388 

TE ( 

6) 

.002997 

.000694 

.000740 

.004166 

TE ( 

7) 

.ool388 

•000231 

,000080 

.002083 

TE C 

8) 

.005774 

.002083 

,001782 

.008333 

TE ( 

9) 

. o o 0 7 «* 9 

.000173 

• 0 G 0 060 

.005250 



F . E . R . 

TE( 1) 
I E C 2 ) 
TE( 3) 
FE( 4) 
TE ( 5) 
T E ( 6) 
T E ( 7) 
TE ( fl) 
TEC 9) 


TflRLE R CUNT. 

COMPARISON Of FOURTH ORDER RETMODS 


KU.-vTZMANN opt. 

4 

k'JTTa 

4 

Sarafyan 

4 

UT 1 

5 

.000277 

.000154 

.000347 

0.000000 

0. 000000 

.001388 

,00?083 

.000001 

0 . 0 no ooo 

,00i398 

. 00 ? 083 

f ▼ 1/ V V JL 

0.000000 

. o o R 3 3 3 

•OOH333 

.008033 

.000006 

. o n 1 1 u 

.000617 

. 001 388 

0.000000 

o.onoooo 

.002777 

.004186 

.000001 

0 . o ft 0 0 0 0 

.001388 

.00?083 

0.000000 

.O0R333 

•00^333 

.008333 

.000006 

. 0 034 r,q 

.002777 

.001041 

0.000000 



TA8LE R CONT 


COMPARISON 


F .E 


:.T ? 
5 

TE ( 

1) 

.000183 

TE ( 

2) 

.001103 

TE ( 

3) 

.000/35 

1 E ( 

4) 

.onl'768 

TE ( 

5) 

.000735 

TE( 

6) 

• 802206 

TE ( 

7) 

. 0 0 7 35 

f F ( 

8) 

. o '' 1 7 6 7 

TE( 

9) 

.000551 


CF FOUHTh ORDER methods 


U T 3 

SHANKS (4-4) 

5 

4 

.000200 

.000006 

.001204 

0.000000 

.000802 

.008125 

.001898 

.008333 

.000803 

.000027 

•002409 

0.000000 

.000803 

,008i?5 

. 0O169R 

.008333 

♦ 000602 

• 0 0 0 020 



TABLE 9 

CCMP ah I SON OF FIFTH OPnF.P METHton<; 



FEHL^FjR 4(S)-l 

FEHLBEHG 4(5) 

f.e.r. 

6 

ft 

T F ( 1) 

.or 1 157 

, 0Ot!4H3 

T E ( 2) 

.on 1157 

• 0 0 id! 4 8 3 

r E ( 3) 

. on 11 57 

.Q02483 

TF ( 4 ) 

. oo625n 

.014743 

T E ( 5) 

.nrlBBl 

,005448 

IE ( 6) 

.001851 

,005448 

TE ( 7) 

.001651 

.005446 

f E < 8) 

.0)6666 

» 0 3 7 1 7 <? 

T F ( 9) 

.001157 

.002483 

T E ( 1 0 ) 

,006944 

.011217 

fE(ll) 

.006944 

• 0 1 12 1 7 

T E ( 12) 

,006944 

.011217 

I E ( 1 3 ) 

.020833 

. q54487 

TE(14) 

.013888 

.022435 

T E ( 1 5 ) 

.01.3888 

•022435 

TE ( 16) 

.013888 

.022435 

T E ( 1 7 ) 

. 04. 1666 

. i 06974 

TE (18) 

. 006944 

• 0U217 

(E(19) 

.001157 

.002483 

If <2'>) 

.00 62 5 0 

, 014743 


FEHLBE»R5<6) 

6 

<?araFyan 

6 

SHANKS ( 5-5 f 
5 

0.000000 

.001111 

.000416 

o.ononoo 

.001111 

.000416 

0.000000 

.005833 

.004156 

0 , ooonoo 

.015000 

.004156 

.00074.0 

.010000 

.013333 

.000740 

.010000 

.013333 

0. 00O000 

.023333 

,033303 

• 0 1 1 n 1 

.076666 

.033333 

0.000000 

.001111 

.000416 

0.000000 

.027777 

.002083 

0,000000 

.027777 

.002083 

0.000000 

.020833 

• 02o R 05 

0 . 0 0 0 0 0 0 

.125000 

. 020835 

,003703 

.083333 

• 066666 

.003703 

.083333 

. 066666 

0 . oonooo 

.083333 

.166611 

.056855 

.416666 

. 166666 

0 . 0 0 0 0 0 0 

.027777 

•002083 

0.000000 

.005B33 

.004156 

0. ocoooo 

.015000 

,004166 



TABLE 9 CCNT . 

COMPARISON OF FIFTH ORDER METHODS 


H'JTOHEP ] 


F . E . R . 

6 

TE ( l) 

. r- n l 1 1 i 

!E< 2) 

. ft o 1 1 1 1 

TE( 3) 

.443333 

T E C A) 

.or 3333 

TE ( 5) 

.006666 

TE< 6 ) 

.or.6666 

r e ( 7) 

.013333 

f E ( 8 ) 

.093333 

TE ( 9) 

. 0 o 111 1 

TE (10) 

.005555 

teuu 

.005555 

TF (12) 

.0)6666 

TE (13) 

.0)6666 

TE(14) 

.033333 

TE ( 15) 

.033333 

TE ( 16) 

.066666 

T E ( 1 7 ) 

. i 66666 

TE ( 18) 

. o 05555 

TE(l9) 

.003333 

T E (20) 

.003333 


butcher 2 



LAWSON 

6 

4 

. 0 0 1 1 1 1 

0 .000000 

.ooiiu 

0 . oonooo 

.003333 

.003645 

.003333 

.002404 

• 006666 

.004146 

.006666 

.004) 66 

.053333 

» 002083 

.033333 

.014583 

.ooiiu 

0,000000 

.005555 

0,000000 

.005555 

0,000000 

.016666 

,018??9 

.016666 

.013020 

.033333 

.02(i833 

.0 33333 

.020833 

•266666 

.010416 

.166666 

.010416 

.005555 

0,000000 

.003333 

.003645 

.003333 

.002604 


NEWTON -COTES Famj, y 


1 uther 

rutcher 

6 

6 

0.000000 

0 • 0 0 0 n 0 0 

0.000000 

0 . 0 0 0 0 0 0 

.008333 

,000260 

,016666 

.001041 

.004166 

.004166 

.004166 

.004166 

.008333 

,002083 

.004166 

.004166 

0.000000 

0.000000 

.010416 

0.000000 

.010416 

0.000000 

.020833 

.001302 

.010416 

.005208 

.041666 

.020833 

.041666 

.020833 

.333333 

.010416 

.541666 

,020833 

.010416 

0,000000 

.008333 

. 000260 

.016666 

.001041 



TABLE 9 CCMT. 

CCmfakisok uf FIFTH dpdEP METHOD*: 



MVSTPOM 

SHANKS (6-6) 

f • E • W • 

6 

6 

T E ( 1) 

. on???? 

.000222 

TF( 2) 

.0022?? 

•0002?? 

TE ( 3) 

.nn 1111 

.000083 

TE ( 4) 

.006666 

0.000000 

TE( 5) 

.006666 

0,000000 

TEC 6) 

. 0 06666 

0.000000 

TE ( 7) 

0.OO0000 

•000222 

TEC 8) 

. 033333 

0.000000 

TE ( 9) 

.00222? 

.00022? 

T E C 1 o ) 

.016666 

. 0 0 1 1 u 

TE ( 1 1 ) 

.016666 

. 0 0 1 1 1 1 

TEC 12) 

0 . o o 0 0 0 0 

.000416 

TE ( 13) 

.093333 

0 .oooooo 

TE ( 14) 

.033333 

0.000000 

T E ( 15 ) 

.093333 

0.000000 

7 E ( 1 6 ) 

0.OO0000 

• ooiin 

1 E ( 1 7 ) 

. 1 66666 

0 . o 0 0 0 o 0 

T E ( 1 8 ) 

. 1 66666 

.001111 

r E c 19) 

. 0 o 1 1 11. 

•000083 

T E (20) 

, o >'6666 

0.000000 


UT 1 

UT 2 

UT 3 

6 

6 

6 

003907 

.000624 

.001244 

0Q3309 

.000624 

.001244 

003908 

.000624 

.001244 

016*08 

.007175 

.002559 

006*44 

.005665 

.001*99 

00**44 

.005664 

.001699 

006644 

.005664 

.001699 

033*73 

.024958 

.006786 

0 0 331 0 

.000624 

.001245 

016412 

.007907 

,008195 

016412 

.007907 

,006195 

016412 

.007907 

.008195 

08?n78 

.009356 

.008622 

03?839 

.015015 

.016393 

032940 

.015815 

.016393 

032R3S 

.015815 

.(11 4393 

164173 

.018713 

.036957 

016411 

.007906 

,008195 

003909 

.000624 

.006226 

014609 

.007175 

,002559 
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TABLE 10 

LOypARI?nN (!F SIXTH CHf'FH METHODS 


F . E . R . 

FEHLRESG 5(6) 
8 

FEH|_HE.-(r, 6(7) 
B 

BljTCWFR (6-7) 
7 

T E < 1) 

. 0 004 op 

0.000000 

0.000000 

IF ( ?) 

.00040? 

0.000000 

.070820 

T F ( n 

. ono47? 

0.000000 

.070820 

IE ( 4 ) 

.or>1587 

0.000000 

,150000 

T F ( 5 ) 

.00 01(15 

0 .000000 

.069400 

T E ( b) 

. 0001 05 

0 . 0 0 0 0 0 0 

,2651 89 

TE< 7) 

.000634 

0.000000 

.265)89 

T E ( 8 ) 

.011746 

0.000000 

3.324034 

IE ( 9) 

. O00402 

0.000000 

1 .133327 

TE (10) 

.on 1005 

.000450 

,054035 

IEU1) 

.001005 

,000450 

,444664 

TE(1?) 

.001190 

.000450 

.444664 

TE (13) 

.003968 

.000450 

3.616)54 

r f < 1 4 ) 

.001587 

.011033 

,397)57 

TE ( 15) 

.001587 

.011033 

2.936067 

T E ( 1 6 ) 

,001587 

.009095 

2.936067 

TE < 17) 

.001587 

.003279 

3.111607 

TE ( 18) 

.001005 

.000450 

6.304039 

TE (19) 

.00047? 

0.000000 

1.133727 

TE (20) 

.001587 

0.000000 

.841 666 

T E ( 2 1 ) 

.00241? 

0.000000 

.04197? 

T E ( 22 ) 

. nn?4i? 

0.000000 

.466975 

IE (23) 

. 00?857 

0 .000000 

.466975 

TE (24) 

.009523 

0.000000 

3.637962 

TE (25) 

.000634 

0 .000000 

.413732 

( E (2b) 

.000634 

0.000000 

2.956993 

^E (27) 

, o o 38 39 

0.000000 

2 ,956993 

TE (2R) 

.070476 

0.000000 

3,159925 

f E ( 29 ) 

.oo?4 1 ? 

0 ,000000 

6.841 975 

TE (30) 

.006031 

.002705 

.327738 

T E ( 3 1 ) 

, 006031 

.002705 

2.237874 

IE (32) 

. Oo7 1 . 4? 

.002705 

2.237874 

TE (33) 

. n?38 09 

.002705 

2.497519 

TE (34) 

. 0095?3 

.066203 

1.931866 

TE (35) 

. o o 9 5 ? 3 

.0^6203 

2,041005 

TE (36 ) 

. or 95 ?3 

.054571 

2,041005 

IF (37) 

.. o ^9S?3 

.019677 

2.622023 

TE (38) 

. oo603l 

.002705 

30.83)620 

T E ( 3 9 ) 

. oo?857 

0 .000000 

6,841975 

T F ( 4 0 ) 

. on9523 

0 ,000000 

9.587962 

T F (41' 

.000105 

C. 000000 

1.027839 

TE (42) 

.000105 

0.000000 

1 . 2 09634 

T F ( 4 3 ) 

.000634 

0.000000 

1.2 0963^ 

T E ( 44 ) 

.011746 

0.000000 

4.83E945 

TE (45) 

.00040? 

0,000000 

l7,07oa?o 

TE (46) 

. OO0809 

0.000000 

1.133327 

TE (47) 

. 005809 

0.000000 

.841666 

TE (48) 

.005809 

0.000000 

.841666 

TE (49) 

.06914? 

0.000000 

235.207200 
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F • F • h’ • 


TF( 1) 
T F ( 2 ) 
T F ( 3) 
T F ( 4 ) 
T E ( 5) 
t F ( b) 
f E ( 7) 
1 E ( 8) 
f F ( 9) 
1 f ( 1 n ) 

f E ( 1 1 ) 
TF (1?) 
T F < 1 3 ) 
TE (14) 
T F ( 1 5 ) 

I F ( 1 6 ) 
T F (17) 
TF(18) 

T F ( 1 9 ) 
IF (20) 

1 F ( 2 1 ) 

1 F (22) 

T F (23) 
TF (24) 
TF (25) 
IF (26) 

T F ( 2 7 ) 

I F ( 2 * ) 
IF (29) 

T F ( 3 0 ) 

T F ( 3 1 ) 
TF (32) 

I F ( 33) 

T E ( 34 ) 
TF (35) 
TF (36) 
TE (37) 
TF (38) 

T E ( 3 9 ) 

T F (40) 
TF (41 ) 

T E ( 4 2 ) 
TE (43) 
TE (44) 
TE (45) 
TF (46) 
TE (47) 
TE (48) 
IF (49) 


TAfel.E lr COM. 


COnP/SPISON OF SIxTh CRpE* METHODS 


BUTCHER (6-8) 
8 

. ''■00661 
.000661 
.001587 
. 004365 
.000264 
. 000264 
.00 158 7 
.00714? 
.000661 
.001653 
.001653 
,003968 
,01091? 
.00529] 
.005291 
. 023809 
.079365 
.001653 
. oo]587 
.004365 
.oo3968 
. 003968 
. o o95;?3 
. 026190 
.001687 
,ool5 8 7 
. o o 9 5 ? 3 
.042857 
. 003966 
. 009920 
.009920 
. 023809 
.065476 
.031746 
.091 746 
. 1 4?857 
.476190 
.009920 
.009523 
.026190 
.000264 
. 00Q264 
.001587 
.007142 
.000661 
.000432 
. 006493 
.006493 
.039069 


SHAnkS ( 6-6)* 
6 

.000619 

.000639 

.000471 

.000698 

.003746 

.003746 

.009783 

.007301 

.000637 

.000636 

.000649 

.000584 

.nnl2o9 

.012624 

.012697 

.022432 

.023823 

.000649 
•000469 
.000698 
.002501 
.002501 
. 0 0333o 

.004190 

.022478 

.022478 

.042031 

, 043809 

.002501 

.002857 

.002857 

.006892 

.007142 

.076197 

.076190 

.141190 

.142857 

.002857 

.003330 

• 00*190 
. 0 0 3 74 6 
.003746 

• 006783 
.007301 

• 000639 
.010800 
.013190 
•013190 
.016190 


SHANKS ( 7- 7) 
7 

.000033 

.000833 

.000489 

,000208 

.002155 

,002155 

.003883 

.003068 

.000033 

.000497 

.000496 

.000022 

.000296 

.000661 

.000661 

.001818 

.001557 

.000496 

.000459 

.000208 

.000196 

.000198 

.001092 
.000203 
.012935 
.012930 
.023299 
.023009 
.0001 98 
.002975 
,002976 
.008184 
.007008 
.003968 
,003968 
.010912 
.009344 
,002976 
.001 092 
.000203 
.002155 
.0021 55 
.003083 
.003968 
.000033 
.000862 
.000602 
.000602 
.000281 



T AbLF 11 


COMPapIsOm CF seventh OWDER METFons 


f.e 


• r? • 


FEWLHERG 6(7) 
1 0 


FEHLREhR 7(d) 
11 


SHANKS ( 7-7)* 
7 


SHANKS ( 7-T ) 
9 


TE ( 1) 
T E ( ?) 
T E ( 3) 
FE( 4) 
TE( 5) 
TEC 6 ) 
TE C 7) 
1 E ( R) 


♦000025 
.000450 
.000077 
O.OnOOOO 
.OOOIBO 
."03156 
.000541 
0 . "00000 


0.000000 
0.000000 
.000066 
•0005S1 
0 .0O0000 
0.000000 
.000462 
• 003«RR 


.0001 15 
.002515 
.001391 
. 0 04*29 
.000578 
.0141 3P 

.oiuil 

.011574 


0.000000 

.000055 

.000183 

.000110 

0.000000 

.000385 

.001286 

.000771 



F • F • M . 

7 E t 1 ) 
T F ( ?) 
T E < 3) 
r E ( 4) 
TE( 5) 
TE ( 6) 
T E ( 7) 
TE( H) 
r e ( 9) 
TE ( 10) 
TE(ll) 
TE<12) 
TE ( 13) 
TE<14) 
IE (15) 
TE (16) 


table i? 


COMPARISON OF EIGHTH ORDER METHOD^ 


FEHLPEPG 7(8) 
13 

.000025 
. nnOD47 
."00088 
.000211 
."00090 
.000166 
."00077 
.001189 
.000205 
.000379 
.000705 
.001690 
. 000720 
.ool3?8 
. 000617 
. 0"95 16 


FEHLBEHG 8(9) 

15 

0.000000 
•000001 
• 0 0000 1 
•000001 
0.000000 
.000001 
. 00000 ? 
.000539 

. oooooi 
• 000002 
.000002 
. 00001 ? 
.000007 
.000006 
.000018 
.00^406 


SHANKS (8-10)* 

SHANKS (8-1?) 

CURT TS 

in 

12 

1 1 

.000025 

.000025 

.000027 

.00149? 

.000029 

.000344 

.001620 

.000029 

.000057 

.00049? 

.000289 

.001762 

.000154 

.000090 

.000100 

.003182 

.000102 

.000079 

• 000?o8 

.000308 

.000196 

.005208 

.002642 

. 0 0 0 8 0 6 

.000205 

.000205 

.000228 

.005327 

.000235 

.002748 

.012962 

.000264 

,000451 

.003943 

.002314 

,014116 

.001234 

.000720 

.000790 

.017319 

.000823 

.000648 

.000247 

. 002469 

,001594 

. 04 1 664 

.021141 

.006500 



CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 

5.1 Conclusions from the Comparison of the Truncation Error 
Coefficient 

In Chapter IV it was stated that the "optimum" set of coefficients 
for a Runge-Kutta method of a particular order could be established by 
determining the set of Runge-Kutta coefficients with the lowest values for 
the truncation error coefficients. Although this is the dominant factor 
in comparing various sets of Runge-Kutta coefficients, the reader should 
also bear in mind the number of function evaluations required by each set 
of coefficients. It was shown in Table 6 that the Fehlberg coefficients 
require more function evaluations for the mth order method when the auto- 
matic step-size control is used. The disadvantage of more function evalu- 
ations is usually offset by an efficient choice of step-size, which in turn 
reduces the number of steps required. Several of the Shanks and Butcher 
methods minimize the number of function evaluations but do not incorporate 
a step-size control. For a rapidly varying function, a poor choice of 
step-size with these methods could result in many more steps being taken 
than necessary, and the resulting total number of function evaluations 
being correspondingly large. 

Since the system of differential equations to be solved determines 
the complexity of the function evaluations and the step-size that can be 
taken, the matter of choosing the best set of Runge-Kutta coefficients be- 
comes problem dependent. For this reason no absolute optimum set of co- 
efficients can be given. If, however, certain classes of differential 
equations are considered, some conclusions can be drawn from the comparison 
made in Chapter IV. 
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Three types of systems of differential equations will be con- 
sidered: (1) systems with rapidly varying functions; (2) systems with 

slowly varying functions; and (3) systems requiring complicated function 
evaluations. For systems of differential equations with rapidly varying 
functions , it is desirable and necessary to incorporate some type of 
variable step control. Efficient determination and correction of the 
step-size is essential if the method is to progress efficiently to a 
solution. The sets of Runge-Kutta coefficients recommended for each order 
method for a system of differential equations with rapidly varying func- 
tions are given in Table 13. If, on the other hand, the system of dif- 
ferential equations has only slowly varying functions, an automatic step- 
size control is not a necessity. The recommendations for problems of this 
type are given in Table 14. Finally, if the system of differential equa- 
tions is such that the function evaluations become quite complicated, a 
method which minimizes the number of function evaluations while still 
producing small truncation errors is required. Table 15 gives the sets 
of coefficients recommended for this type of problem. 



Table 13. 


Recommendations for Rapidly Varying Functions 


Order 

Recommended Method 

No. of Function Evaluations 
Reqd. with Step-Size Control 

3 

Fehlberg 3(4) 

5 

4 

UT 1 

6 

5 

Fehlberg 5(6) 

8 

6 

Fehlberg 6(7) 

10 

7 

Fehlberg 6(7) 

10 

8 

Fehlberg 8(9) 

17 


Table 14. 

Recommendations for Slowly Varying Functions 


Order 

Recommended Method 

Number of Function Evaluations 
Required 

3 

Fehlberg 3(4) 

4 

4 

UT 1 

5 

5 

Fehlberg 5(6) 

6 

6 

Shanks (7-7) 

! 7 

7 

Fehlberg 6(7) 

10 

8 

Fehlberg 7(8) 

i 

13 
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Tab.le If). 

Recommendations for Complicated Function Evaluations 


Order 

Recommended Method 

Number of Function Evaluations 
Required 

3 

Ralston Optimum 

3 

4 

Kuntzmann Optimum 

4 

5 

Shanks (5-5) 

5 

6 

Shanks (6-6) 

6 

7 

Shanks (7-9) 

9 


8 


Curtis 


11 


5 . 2 Topics for Future Study 

With the availability of the equations of condition, it is pos- 
sible to attempt to produce explicit Runge-Kutta methods of the ninth, 
tenth, and even possibly higher orders. However, the solution of the equa- 
tions of condition for these high order methods will be very difficult. 

It is hoped that a computer process can be developed to solve these ex- 
tremely large systems of nonlinear algebraic equations. Aside from 
producing new methods, the equations of condition can be used to improve 
the lower order methods. The UT 1 coefficients were successfully developed 
previously to optimize the fourth order Runge-Kutta method. If a computer 
process is developed to solve these equations of condition, then some form 
of optimizing process can be used to produce similar sets of coefficients 
for the higher order methods. 



APPENDICES 



APPENDIX 


Original SYMBAL Program 


SYMBAL VERSION 1.1E 1 3/14/69) 

START PROGRAM ", 0, 0 

equations of condition for runge kutta mEthoo 
=begin= 

rNEWE DX»DY*Ym .YN?*TS t ZZ,CCEF| 

Du^P Is 51 
H *■ 41 
P Is 3| 

m 0DEC4] is 

F i * i*OI = F0R= I t ■ 0 1 H- 1 =CO= »0 im-i-1i++| 

T5 is EFORE I Is 0 | m~ 1 ESUME EfOhE j i* Oim-1-j ESUME 
H*OXtI*DYtJ*F[I f J]/FACf (I)/FACT(J) | 

A * s *P*4| 

K I ■ *0 lH*Ft J tPl4* 

tile* = FOR = I is 1 j P EDO= f* 0 | T - 1 ; * * { 

=FOR= I ** IIP =D0= 5REGINE 

S !■ EF0R5 J «■ Oil-l ESUWi BlI,jHnaj]l 
Mil l« =s= 0 X IA[ I]*H,O y is S TS* 

EENOE » 

C 1 s *o * P * * I 

YnI 1* EFORi I Is °IP ESUKr Cl 1 1*K 1 1 1 ’ 

2 I c sF[0 t 0) t M|*| 

YN2 is H # Ztl]*( EFOR= I l« 2im ESumE r «* \/f ACT ( t ) • (Z t T J is 
EFORE J I* 01 1-2 =SUH= EFCR = N is 0 1 1-J-2 ESyM = 

( Ff J4l tNj^Ff Of 0 J«Ft JiNMJl* ED= F C J S ZtT-i)>>» 

11 is SPLIT POWERS (Yn2-yn1*H)» 

EFOR= I Is 1 1 M =DC= EFORE J Is 0 1 1-1 EDO= EFOR= N is OlI-J-1 =00= 
ZZUl Is E§= Ff J,N J iQtFt JtNl % Z7 C I 3 1 
EFORE I Is IIm EOOE EBEGINE 

t is sp l it powers <zzt n t c> » 

EforE J 1 s 1 SIIPBOUNO (T) =COE =BEGIN= 

CnEf IS TC-JH 

t cjj i* T c jn 
Eenoei 

EENnE * 

SEnD=I 

012277R 
0. 0.642 



total space Available 
execute °RogRam 


l 


I 

Revision 1 


CV». )tl A(_ VFQ l,|r I.T/l 

SIAM ptioncAi.- n, r. o 

T hf FCti*TICN«; r >r Cnwnn . F(> v GF» rRAu" 'cij* Gf-ki itT a -icTh'D 
Foo *- N: Y '>■ F = NL» G" *"l T1 A f Fo* A I iTONr w m iC SYqTF • 

=fEMn = • 

= kf;.,= n*. r Vf Y N 1 » y.rt T . // . f " F F ; * 

p : = 4 i • 

P : = U 

T AY LOP CFRIF.S M CftEf4i •- ' j . 

F I s *n: y-l I * ; 

CO*FUT,ttOK CF W T.-Vi nR~SF»IES FOl ? Tnr 1 k : F ^ ™ ’ ‘ 

B t* +*F OP z T ; * 150 =uO= *011-1:4*1 

EFOR= Tl* HP EUnE EBEGINE 

c t» E Fr R = j.s n i t — 1 Ecl.m = R r t , j ] *k f j -j : 
will i* ES= OYiS $ tc ; 

EFN.-ri 

FINAL RiinGE k \ T T'/* S»|M»VTI0N| * 

C : * -a : pi*i 

Y N 1 := :FCP= I Is 0 1 P EVJ*E C[Il*srHI 

7 I » -»f f M : + I - • 

* K? l^r #ZCn 4 (SfpnE ; is 2 1 m Eci:m= Ht T /FArT(T)*(Zrn is * 

~ J,a ,J5 I*2 r?U ■*= < F c 0 ] «F r IM J) «EP = Ff Ji « 7rr-llni 
CIFFFuF.,rE JF yni a n u v^2 SFT To 7Ep« YTtLr<; rora T r^i<; : 

27 Is SPLIT POWERS { ) t 

ifra l zzT, J; M HD ° 5 - rCH; * is 0!t -’ = u0; 77tT1 •' frJlTO.rt-JJ ! 

SEcAPSTTp, r;F EquA T TQ‘ S ^CCocCIvG Tq pfloUUCTS n p n 7FrF-fP N T F si | 

Di"’P ;= c. 

r F - k= t := i E dge epeginE 

T i« SPLIT Pon!-p5 (ZZI T 1 *0) I 4 

=PO«E J la IjUPTOLNO(T) EDGE =«FGIN= 

CCFF :* Tt-jli " 

TC 1 := Tt in 
=FNCE; 

=ENn=i 

SfNr=! • 

• *#«>*•***<.«»******, »«««•«» t » 1 

TOTAL SPACF fivATLAPLF p2i 277« 
r*EC».»TE PROGRAM f . . ^,0^7 



-P 

O 



Revision 2 


SYMdAt VERSION 1 . 1 F i 3 / 1 4 / 69 ) 

STAHT program n • A. 0 


THE FCIJATIOkS of CONDITION FOR A GENERAL RUKGE-KUTTA mfThOO 
=BEGIN= 

rNFWr OX* OY* Y*M * YN2 * T5* ZZ» COEF l 

h:*s» Hl*5l 

MODF{4) Is ‘'-I * 

M0^ F t5] Is 01 moOEM 3 l* Ml 

F is 

* D I * •* 1 IH-1 1 M 

G IS r* 1 1 M-l 1 + 1 

► T$ 1« EFOR= I :* i:m- 1 =SLM= n«|) [ M *F ( IJ /f ACT { T ) I 

Ts i« H«F( 0 l*TSJ 

Computation of t*f taylor series foh the ki 

K Is f»OtH«Ft0J »Rt At 

B Is •*=FQR= Its UP =DO= HOll-ll**! 
sFORE II* UP =DOr =BEGIN= 

MODEU) is M-21 

GC i 3 I • SFOR= J l« otl-l SSUM= HtI,jJ*KUll 
=FO«= N la 2 : m-i =DO= =8EGlN= MODE [ * ) is m -N- i | 

OCN1 • * ocn*i t 

MODE ( 4 ] is M.2I 
GtNl I ■ 0tN)+N \ 

=FNO= I 

MODE 14) I* m-1 I 

K(U Is =S= DClltGlil , 0 1 ? J : G 1 2 J * 0{J)lG[3l * 

0U1IGC4) * TS I 

“Eno=* 

final PUNGE KUTTA summation: 

C I * *0lPt + l 

YN1 is =FOP= I 1 * OlP =SUH= cl I III 
Z ts..*FtO]«Pj*l 

YN2 *s m*Z 1 1 ) ♦ (=FOR= I *« 2 1 M =SUM = H4J/FACTU)*(ZU3 

=FOR= JIb 0 1 1-2 =SUM= <Ft03 # FC J*1 1 >*=0 = Fm % zti-Dm 

sENOs! 


TOTAL SPACE AyAlLARL^ 0\2*T?fi 
EXECUTE PROGRAM 0. 0.6U 



Shanks * Method with 
Fifth Order Output 


SYMbAL VERSION 1 * IE 
STftRT PROGRAM 


=begin= 

DUMP 1*0 I . 

MOqeMJ i* o » 

Ml.51 

NlaBI 

Lla4^|l|lt2 f 4 t 9t20,49fll7»297»74#»,10i4* 

A|,itl|N|*» 

C f 1 INUI 

Bl *#5FORE Its 1 1 N EDOr *0*1-11**) 

Ql «**FqR= K !■ llM-1 E°OE 
#SFOR= I j ■ 1 l N SOO= 

♦1 tLCKl 14**1 

COMPUTE THE Values Of THE C s I 

= F0Rr I lal IN EDO? Q C 1 * I * 1 3 la A t 1 3 » 

=EoR 5 KI.2IM-] =0o= 

=REGlN= 

=FOR= 1 1 a 1 IN =00= =REGIN = 

=FOR= TlclllCK-1) EDO; EBEGINE 
Q[K*I«T3 l« A[l] # OtK-l*I»T)* 

OfK*I»T*L[K-l3) la K* (EFORE Jlslll-l 
ESUME (B[ I. J3*Q[K-1«J|T3) )* 

SEND- I 

r IFE K <4 5THEN= =GOTO= CUIT | 

Pl«l-C*-3M u:*2*L£K-l ) I 

rFORE T 1=^1 : l, *L[K-2 J-LCK-3] =00= =dEGlN= 

R I *R* 1 I 

OCK.ltT) la2*(=F0R= JlaltI-1 =SUM= B[ItJ)«A[ji> 
♦QlK-2*I*Rll 

5EN0=» 

=IfE K<6 ETHENE EgOTOE quit* 

=FoR= vial 12 E°0= =8EGin= 

R> “L C K-4 1 $ Ul*T-l» 

5F0R= T I *U* 1 I U*LtK-3)*LlK-4)-LtK-5)*llK-6l 

5D0E EBEGINE 

R)»R*1 J 

OlK*ltT3 Ia3* (=FOP= =SlJM = 

BtT,Jl«0[2,J,V))*Q[K-3 f I*Rj| 

5END= I 5e>*d = I 
5lfH K<8 ETHENE EgOTOE QUIT! 

-FOR? VJ« 1 1 4 =00= EBEGINE . 

Ri*LCK-5Jl Ul*T-i» 

SFORE T I » U*1 I U*UK-A}-L[K*5l-l[K-6)aLtK-7J 

=oo= =begin= 

RloRtll 

0 [ K * 1 1 T 3 la4a( = F0R= J I a 1 1 1 - 1 =SUM= 
B[T f J]#Q[4 f J,V])oQ[ K.4, I*R) | 

SENDS f SENDE ) 

•CUIT* 

SENDS * 

. SENDS * 

2F0R= K *• 1 IM-1 =00= SBEOINE 
xFOR= T Im 1 1 L £ K 3 EDOE EBEGINE 
Dump u o i 

SBEGINs 

X Is EFORS I la 1»N =SUM= C t 1)* < 01 K 1 1 * T ) ) » 


0 


PlJMP | a5 | 

0 : * 1 / < k ♦ 1 ) - x i 
OU M P (=01 
EEMPE i 
=FNP= * 

S E N n = j 
=ENP=| 


total space available o?i3is« 

execute PROGRAM 0. 0.955 

Ct2]*A(2] - Ct33*At33 - Ct4)*A£4J 


1/2 . c r 1 1 * A C 1 3 - CI2]*A(2j - c£ 3 )*a£3] - c t A J «A t A ] 

- c t51* A t51 I 

. P# 2. 2R 

l/ ? 7 " C[23*At?Jt? - Cl3j«Af13^2 - C ( A ) a A t 4 j *2 

- C l S | 

„ 0. 2. 75 

1/3 - 2#P[2,n*C[?JaA[li - 2«C[ J j*(8( 1» 1 1*A[ 1 ] 

♦ 0C3*21*Af21> - 2*C[Al*<R[4,n*Atn ♦ Rr4.2l*Af23 

♦ bt4»3]*A[3j> * 2*C(5 ]« (R(5* l ] aA [ 1 3 ♦ P[S*?1*A[2) 

♦ b[ 5* 3 ) ®A r 3 3 ♦ Rt 5*4 J*A( 4 ) ) I 

0* ? • 1 2R 

1/4 - C( 1 )*A [ 1)43 - C[2)*A(2i+3 - C[3]«Ar3)+3 - C[4)*A[4)*3 

- Ct5J«Af5)A3l 

*. 2,174 

1/4 - 2*B[2*1)*C[?)«A[1 ]*AI23 - 2*c t 3 1 * a n 1 * <B £ 3 * 1 1 « a Cl 3 

♦ B[3*2) •*(?)) - 2*C{4 ]*a [ a j # { e [ 4* J ]*A [ j ] , H[4.23*A[?3 

♦ 0 £ 4 • 3 ] #A( 3 ) ) - 2*C(5j*A[5]#{B{b,l3«A[ i) 4 B[5,21*A[?3 

♦ bt5*3)*At3) ♦ RC5.4 )*A(4 J > » 

1/A - 3*Rt2»lJ # C(2]*AClJt2 - 3*C{3)*(Bf3.i )*Atl 3*2 °* 7 ' 22 * 

♦ b 1 3 #2 J*A f 2 3 12) • 3«C C4 1* tb (4 ♦ 1 J«A [ i J +2 ♦ r [4 » p JaA C? 3 *2 

♦ ot A .3 ]«A{ 3)42) - 3«C[ 5 ]* {H[ s, 1 3*A[ 1 3^2 ♦ P£ 5 « 2 1 • a [ ? i ♦a 

♦ bt5*3)*A[3l42 ♦ BC5fA )*A£4142U 

1/A - 6*P(3t2]*Bf2,llaCc3J*Am - 3*C C 4 ] « ( ?#B { 4 , 2 3 *B r 2 f fj ^ ? 
*aU1 ♦ ? *rC 4 ,3)*(r{3,i 3 • 4 1 1 1 ♦ bf3*2 J«Ar;>l) ) 

- 3*C(5]* ( 2*B[5 # 23aH(2*l ]*A[ t j * 2*B[ 5, 3 ] * ( R[ 3, ^ 3#A(l 1 

♦ bt3*23*A[2)j ♦ 2*Bl 5,43# <0 [ 4 , 1 J«A[ i J * Pr4,2]«Ac23 

♦ b[4*3]*A[ 3) ) ) 1 

<- „ - ^ P* 2*343 

I/ 5 - C[l)*AfiJt4 - C(23oA[?^ + 4 . C[33«Ari}44 . C[4]*A{4]t4 

- Ct5)*A[^)*4| 

„ n. ?,39i 

)/f * 2*B[2.1 3#Ct?)»A£ 1 jaA[2)42 - 2«C c 3 ] • A f 3 ] 42* ( B [ 3 , 1 3 *A [ 1] 

♦ d ^3 12 )*A 1 2 3 1 - 2*CC4lnC4)A ? «(riC4 M )«A(i1 ♦ rU* 2 J*a£ 2] 

♦ btA f 3]#A{3J) - 2*C£ = 3*A [5l42* (0(5 f j j*a n 3 ♦ H[5,2]#a[2] 

♦ HC5t3l*A[3J ♦ Bt5,41«At4J> I 

_ 0* 2*487 

1/5 - 3*Bt2*l3*CC2]«ACnt2*A[21 - 3*C [ 3 1 *A { 3 ]# (Bt 3 . 1 1*A[ 1 ) 42 

♦ bE3*2)*A[2)4 2 ) - 3*Ct43*At4l*(dt4,i JoATjla? ♦ « £ 4 , 2 1* A 1 2 1 ♦? 

♦ 0t4t3]»A[3]42) - 3»Ct5l*AC51* (b[5, UaACi ]42 ♦ B£ 5 , 2 i *A r 2 3 4? 

♦ «C5*31*A[3l42 ♦ B£5*4 J*At4l42) ; 

, P* 2*547 

1/5 * 6«P£3.?3*H£2f 13 «C [ 3 ]• A [ 1) « A [ 3 J - 3«C f 4 3 * A [ 4 ] a ( ?*B [ 4 , 2 } 
*Ri2,i1*Afjl ♦ 2*Bf4t3l*(R.t3,i )* A ti 1 ♦ 0 £ 3 , ? 1 *a t 2 ]) j 

- 3*C[5]aA C 53o ( 2*B[5,23*P[2 f 1 3*A[ 1 3 ♦ 2*P £ 5 , 3 ] a ( B[ 3 , 11 a A [ 1 3 

♦ b£3.23«A[23) ♦ 2aHf5*43«(0(;4, 1 j*A[ \ ] * Rr4 1 2]aA[2] 

♦ d£4*3)*A£3] > > I 

c n. 2*605 

1/5 - 4at-[?* ] ]*C[2]«A[ 1)43 - 4aC [ 3 ] a ( R( 3 , ^ -j »a [ 1 ] 4 J 


•P 

ro 



Shanks’ Revision 1 


o 


o 


0 


0 


♦ Bt 3»2 )*A r 21 *3) - 4 «CC 4 j#(Hr 4 » 1 1 1*3 ♦ R I 4 , ? J *A[ 2 1 *3 

♦ BU*31*M3] + 3> " 4*Ct5l*<H[5*l )**l\U3 ♦ R[ 5* 2 1 *A( 2 3 13 

♦ bl5*31«A[33t3 ♦ B15*4)«A14]43>| 

6 * 2*665 

I* 1/5 - 0*H[3i 2]*HI2»1 I«C[3J°A{ l ]*A(2j - 4*r [ 4 3 • ( 2*b [ 4 , 2 ] 

*B&2i 1 1*A h 3 «a t 2^ ♦ 2«8 ^4 *3 J *A t 3 )• (Bt3 t i l*a T i 1 

♦ Bt 3»2 1 * a 1 2 1 > ) - 4*c[5 ) # ( 2«RC 5 *^1*B(2 t l l*Af n«A[?l 

♦ 2*BC5*3J*AC3]4(B(3«1J*A[1] * B[3,2l«At2l> 

♦ 2*B[5,41«AC4]<MB[4,l]«Atl] ♦ B[4.2]*A[2 i ♦ R [ 4 *3 )«A (3 ] ) ) I 

0, 2,726 

Is 1/5 - 124BC3 f 2)*BC2 # U«a3)*AC 1J + 2 - 4©c r * 1 • l 3»B [ 2l»Bt *t 1 1 

♦ A113*2 ♦ 3*Bt4*3l4(BC3*l)4A[i Jt2 ♦ 8 C 3 ,2 1©A [ 2h2 > ) 

- 4 # CC5]*(3 # Bt5i2)*8C2*l]*A[i)t2 ♦ 3©B C5 .3 1 * <R 1 3#1 ) »A j 1 1+2 

♦ BI3*2 J«A{23*2> ♦ 3«BC5,4J«(fl[4»l)*AClU? * Bt4.2)»Ar2)*2 

♦ B[4* 3]*Af 3)t2 > ) \ 

6 * 2*786 

l> 1/5 - 24 «Rr 4 , 3 ]*B[ 3 ' 2 )*e[ 2 «l ]«C[43*AC 1 ) 

- 4 # CC51*<6*BC5*31*BC3*2]*BC2*ll # A(il ♦ 3*n f 5 * 4 3 * < 2*R f 4 , 2 J 

•Bt2,ll*Ati3 ♦ 2* e t4,3)« (Bt3, 1 1« A C 1) ♦ R M . ? 1 * A [ 2 1 >) ) l 

ft. 2.833 

:> 1/5 - 4*C£23*tb(2,l l«A(l3>t? - 4*CC 3 !• (0 fl. U»AC1 1 

♦ B t3,2 J#A(2)> t2 - 44CUJ©(BC4,l3*A[j] ♦ Bt4.2l«At2) 

♦ «14,33©A(33)42 - 4©Cc53#(B[5,i3*A[i) ♦ «r5.23*At2) 

♦ Bl5'33*AC 33 ♦ Bl5»4J«At4))42l 

ft* 2*893 

EXECUTION COMPLETE ft. 2.897 


SYHOAL vt^Sln* 1 • 1 r 
STAHT PPOGRaw 

EbEGlMr 
OUmP : a 0 I 
MODE ( 1 J »• 0 ; 

*• ; s h 1 N r 3 1 2 X 

Ll=i»0tltltl«l«1.2*4fH.l6t3i: + i 

A! 3 . j |i«L| 

c 1 o * 1 : 1 * 1 

Bl»r*5FJH= I «s 1 1 N =00= *0ll. ll*4| 

EFORE =DOE RlI*lM a 0» 

EFoR = =Doi Hll* 2 ]:aft| 

EFORE 1 1 *5 1 n = 00= R{ l * 3 1 : »0 l 

0(«4 = FQR= K : « 1 : M-l =00= 
f* =F0W= 11 b UN =00 = 
f*llL[K] 1 * 4*1 

= FoR= 1 * ■ 1 * N =D0r OUfI.ll*«At!l| 

=F ORE IIbUn =00= Qt 2» 1 1 1) I "A [ 1 1 ©2 I 
=F0«= I t«l|N =00= Q( 3 • 1 1 1 3 1 *A [ 1 1 *3 1 
=F0R= 1 1 * 1 1 N =005 QUtI*l)l«AtI»*4l 
5FOR5 K * »5 ih- 1 =00= EBEGINs 
5FOR5 =00= =bEGIN= 

EFORE T«b1il[K- 11 =00= EBEGlNE 
0CK.I*TJ1»A{ I]*U[K-ltI*T J * 

=EN0=» EFORE T : = L i R- i 3 ♦ l I L [ K J =00= 

OtK.I *T3 t«K* ( = F 0R= = Sijm = <0 1 i j ] *C r k-1 * J. T-L Ck 

EEMD5 ♦ 5EN0E 1 
=FOR= =00= =BEGIN= 

=F0R= MsliLtKl E 00 E ErfEGINz 
OumP i—O ( 

EFEGINE 

Xts =F ore i : a 1 : n ESomE CC 1 1 * (Ql*. i .T J I \ 

FILE2 »« EEFUNCHEE ? 

= IF = K = 7 =ThENe OUMP j = 5 i 
Ola 1/(K*1)-X 1 

=EKO=* =END= ? EEuO= * EENCE I 


TOTAL SPACE AVAILABLE 
EXECUTE PROGRAM 


1 3/ U/bQ) 

A , (> , t.) 

o a o * a 


-ill)); 


122315b 
ft. 0.671 


-P 

CO 



APPENDIX II 


All of the coefficients used in this study are given in rational 
fraction form since the SYMBAL language uses only rational fractions for 
computational purposes. All of the coefficients with the exception of 
the Fehlberg 8(9), Curtis, UT 1, UT 2, and UT 3 sets were originally 
developed in rational fraction form. The exceptions were developed in 
decimal form. For this study, these decimal coefficients were taken to 
six decimal places and put over a denominator of one million to give the 
necessary rational fraction form. 
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third order coefficients 


FEHLBERG 6 3(4) 

C [0] ,*79/490, C11J,*0, Ct2J ,*2175/3626, C 13] 1*2166/9065, 

A 111 1*2/7, A 121 l *7/15, AI31l*35/38, A[4),*l, B[l,0J,*2/7, 

B[2,0J, *77/900, B[2, 11 , *343/900, B[3»0J ,*805/1444, B C3 , 1 1 , — 77 1 75/ 
54872, B [3#21 ,*97125/54872, 

CLASSICAL 11 

CI0JI.3/8 , Cini.2/1 , C121H/4 , AtDl-1/2 I t!2)l>t I 
B tl #0] 1*1/2, B[2,0]| — 1 , Bt2,lJ|s2 , 

RALSTON OPTIMUM 15 

C [0] , *2/9 , CI1J ,*1/3, C 1 2 J 1*4/9, Atl]|*l/2, A[21,*3/4, 

Btl#0J 1*1/2, B- [ 2 » 0 J |*0, B[2,l),*3/4, 



HEUN 


7 


C C03 I * 1/a f Ctm®0> C[2]»«3/4y Atni = l/3f At2]**2/3» 

B [1 *01 1*1/3 t B [2,0] 1 = 0 » Bt2*l)tB2/3> 

NYSTROM 14 

C(0]|sl/4| CClJta3/8f C C2J 1*3/8 y Atl]l«2/3f At21t*2/3y 
B(l,0]|«2/3y B (2» 0] I »0 y Bt2,tl|B2/3y 


P 

CD 



FOURTH ORDER COEFFICIENTS 


FEHLBERG 6 3(4) 

C (0] 1*229/1470; CtlJ|«0? C t2J |»1 125/1613? C C31 tail371S/61585f 
C l«l t a 1 / 1 8 1 

A 1 1 ] I °2/7 ? A 12J I »7/15 ? A (33 t ■35/38 f AtdJt«l? Btl,0J»=2/7? 

B t2#0] |«77/900| B [2, 1] 1*3(13/900* B [3 $ 0] I *805/ 1 444 » B (3, 1 J ia-77175/ 
5(1872? B 13,2] 1*97125/5(1872? 

B (4,0] 1*79/(190? Bl(J,i]?»0? B t4,2 J »b2175/3626? 

B U#3] 1*2166/9065 ? 



FEHLBERG 5 fl(5)-l 


C10] I * 1/9 | C til I® 0 f C121 »« 9/20 f C [3] |® 16/45 » 

C[4] »® 1/12 f 

AtlJ I" 2/9 f A 12J 1 8 1/3 j A [3] I® 3/4 > A [4J |« 1 f 

Btl»0] I® 2/9 | B[2,01 |* 1/12 } Bt2,ll I® 1/4 f B[3,0) I® 69/128 » 

Bt3»l] i® -243/128 f B[3,21 |® 13S/64 ) B[4,0J |« -17/12 f 

814, 1J I® 27/4 i Bt4,2J I® -27/5 } Bt4,31 *■ 16/15 ) 

FEHLBERG 5 4(5)-2 

Ct01 I® 25 /216 | C til I® 0 f C 121 t« 1408 /2565 | 

C131 |® 2197 /4 104 i C 141 »«-l /5 » AtlJ I® 1 /4 » 

A [21 I® 3 /8 » A [31 t®12/13y A [43 t « 1 f Btl,0J I® 1 /4 j 

B[2,0J i® 3 /32 > B [2, 11 1® 9 /32 t B[3,0J |® 1932 /2197 » 

813,11 I® -7200 /2197 i B[3,2] |® 7296 /2197 » Bt4,01 I® 439 /216 | 

8 14,11 l® -8 t B [4, 21 I® 3680 /513 | B[4,31 »® -845 /4104 » 

-P 

CO 



CLASSICAL 11 


Ct0]l = l/6> Ctllisl/Si C [2] I » 1 /3 j C £35 I -1/6 > AtlJial/2, A(21t«l/2 , 

A [3] I * 1 1 B [ 1 , 0J I »i/2 f B 12# 03 I >0 j B [2# 13 1*1/21 B(3#0)|>0j B[3#13i&0f 
B 13# 23 l*lr 


KUNTZMANN OPTIMUM 9 

C I0J |«55/360» C (11 |>125/360> C 123 l>125/360| C (33 t >55/360 f A C 1 1 I «2/5 f 
A (2] 1*3/5 j A £31 I >1 > B(l#03i>2/5| B £2 # 03 t»-3/20? Bt2,i)|«3/4f 
B (3#0J »>19/44f B (3# 13 |*«15/44j 8 (3,21 »>40/44» 

KUTTAiO 

C(03l>i/8» C C 1 3 l»3/8» ct23l>3/8| C[33t«l/6' A(l)t>l/3* AC2]|>2/3? 
AC31|slr B £1 # 0 } }>l/3f BI2#03t»l/3| 0 [ 2 # 1 3 t « 1 f B(3#03t«|f 
B 13# 13 »«•!» B [3# 21 i«l> 



5ARAFYAN18 


C[0J I* t/6 t CCtJ is 0 f C [2] is 2/3 , c 131 is 1/6 , 

All] is 1/2 , A [2J I s 1/2 f A [3] is 1 , 

A t4] is 2/3 , 

Btl,0J is 1/2 , B12,0] is i/a , B12,l] is i/a , B13,0] is 0 , 

B13,i] l« -1 f B (3 , 2] is 2 , 

Bta»0] IS 7/27 , B[4,1J IS 10/27 , B[4,2] |s 0 , B14,3] is 1/27 , 

UT t 1 

C 10] 1*124584/1000000, Cll] is0| C f2J 1*577732/1000000) CISJis 
1 1 167899/10000001 C 14] I ••108702 16/ 1000000 1 

A til 1*266000/1000000, A 12] 1*398999/1000000, A 13] 1*997475/1000000, 

A [4] 1*1, 

B [1*0] 1*266000/1000000, B 12,0] 1*99749/1000000, B 12 , 11 1*299250/ 
1000000, B 13,0) l*997397/1000000» B 13 , 1 ] | *-3740204/ 1 000000 , 8 13,2] |* 
3740282/1000000, B (4,0] 1*1018551/1000000, Bt4,l] I *-3826725/ 1 000000 , 
814,2] 1*3810768/1000000, B14,3] I *-2594/ 1 000000 , 



C 10) 1*125181/1000000) Clll**0) C [2] |O299910/1000000| C [3] I *437933/ 
10000001 C 141 1*136974/1000000) 

A IU I *203066/ 1000000 1 A 121 I *304600/1 000000 1 A [31 I *620351 /\ 000000 j 
A (41 1*1) 

B [1,01 1*203066/1000000) B 12,0] 1*76150/1000000) B (2 , 1 1 t *220450/ 
1000000) B 13,01 1*327625/1000000) B 13, 1] I *« 1 0 1 6945/ 1 000000 ) 

B [3,2] 1*1309671/1000000) B [4,0] I **300307/ 1000000 ) B (4, 11 1*2751166/ 
1000000) B [4,2] **-2664664/1000000) B (4, 3] I * 1 21 3805/ 1 000000 ) 

UT 3 1 

C [01 1*128196/1000000) Clll^*0) C (21 1*275745/1000000) 

C 13] 1*454937/1000000) C [41 **141120/1000000) A (lj 1*200933/1000000) 

A [2] **301400/1000000) A [31 **606170/1000000) A[4}**1) 

Bll,0] **200933/1000000) B [2, 0] **75350/1000000)0(2,11 **226050/1000000 
) B [3,0] 1*308207/1000000) B[3,l] **-934792/1000000) 

B 13,21 **1232755/1000000) B 14 , 01 **-232399/1000000) 

B [4, 1] **2571843/1000000) B [4,2] **-2609052/1000000) 

B [4,3] **1269608/1000000) 



SHANKSU-4) 


C t0] t*«179124/70092r C tlJ I >200000/70092 t C 12J l*40425/70092| 

C t31 1*8791/70092* A(l)l«l/ 100 | A[21|*3/5| A C33 l« t f B tl ,0] |al/l00| 
B [2*0] t*-4278/245f B C2 r 1 1 1*4425/245 f B 13,0] 1*524746/8791 f 
B[3,U ta-S32125/879l» B13,2] 1*16170/8791 » 


FIFTH ORDER COEFFICIENTS 


FEHLBERG 5 4(5)-1 

C 10] »s47/450» C tl] |s0| C (2] l*12/25| C 13] 1*32/225 f C(4]|sl/30f 
C [5] l*6/25f 

All] I* 2/9 f A 12] |s 1/3 } A [3] |« 3/4 » A (A] |« 1 j A 15] t *5/6 9 
Bll#0] |« 2/9 f 0(2,0] is 1/12 p B (2, 1 ] is 1/4 f 0(3,0] !■ 69/128 i 
B(3,1J is -243/128 j B(3,2] ts 135/64 | B(4,0] is -17/12 $ 

B(4,1J ts 27/4 | 8(4,2] ts -27/5 | 8(4,3] is 16/15 » 

B (5,01 »s 65/432| B (5, 1] »»-5/16f B (5,2J I>13/16| 
B (5, 3] l«4/27 1 B (5,4] t«5/144» 



FEHLBERG 5 4(5)-2 


C [01 |sl6/135| Ctljl*0» C t21 I *6656/ 1 2825 » C [3) I «2856 1 /56430 f 
C [4] I ■•9/50| C [51 l «2/55 » 

A [1] |sl/4f A (21 I »3/8 > A[3]|ai2/13f A[4J|sl, A[5J|«l/2f 

B [ 1 » 0J |s 1 /a | 

B[2,0J »b 3 /32 i B[2,1J l« 9 /32 » B[3,0) I* 1932 /2197 f 

B(3,ll l« -7200 /2197 j B[3,21 |« 7296 /2197 » B[fl#0) !■ #39 /216 | 

B [4 # 1 J is .0 f B [4# 2] is 3680 /515 » B[4,3] l= -845 /4104 f 

B[5,0] ls-8/27 1 B(5,ll t«2f 

B [5,2J ts-35aa/2565| B [5,31 »s!859/4104| B [5# 4] »=» 1 1/40 f 

FEHLBERG 6 5(6) 

C[0] l» 31 /384 | CIl) l« 0 f C[2J !■ 1125 /2816 | 

C [31 is 9 /32 | C 14) I s 125 /76B j CIS) is 5 /66 » All] I* 1 /6 » 



At2J IS 4 /15 » Am IS 2 /J p A [4] is a /5 » A [5] |s 1 , 

Btl,0] is i /fc p 0(2,0] is a /75 | Bt2,lJ t s U /75 p Bt3,0]|* 

5 /6 | 013,1] is -8 /3 » Bt3,21 |» 5 /2 p BI4,0] |« -8 /5 p 

Bt«#n is iaa /25 p Bt4,2] is -a p Bta,3i i» u /25 p 

B 1 5 , 0 ] |S 361 /320 p B (5, 1 ] is -is /5 p Bt5,2] |* «07 /128 p 

Bt5 # 3] is -ii / 8 0 , 815,4] is 55 /128 p 

SARAFYAN18 

C 10] |sia/336p C [1] |s0p C t2]ps0j C (3] |&35/336pe C4] t*162/336p 
C 151 I*l25/336p Atl]|si/2p A(2]|sl/2p At31i*lp AJ41|s2/3p A(5]|s2/I0p 
Btl>0]|sl /2p B! 2 , 0 ]psi/«p Bt2,lll*l/4p BC3,01ts0p Bt3,l]l«-lp 
B C3,2] | s2 p B 14,0] |s7/27p B (4, 1] i»10/27 p Bt4,2]|s0p B [4# 3] 1*1/27 p 
B [5,0] l*28/625p B C 5 , 1 ] I s-1 25/625 p B 15, 2] 1 8546/625 p B 15, 41 IS-378/ 

625 p B C5, 31 l»54/625p 



SHANKS15-5) 19 


C [0] *0105/11341 CtlJ»*0» C 121 ,*500/1134, C 13 ] l *448/ 1 1 34 , 

C [4] 1*81/11341 C15]»«0, A 111 1*1/9000, Al2],s3/10> AC3],*3/«, A[4]|*l 
, B [1 ,0] 1 = 1/9000, B [2,01 13-4047/101 B 12, 1] 1 = 4050/10, B[3,01,a 
20241/8, B(3,l] 1**20250/8, B 13 , 21 1*15/8, B 14,0] ,*-931041/81 1 
B [4# 1) I "931500/81 , B 14,2] ,*-490/81 , B 14,3] ,*112/81 , 

BUTCHER 2 1 

C10J 1*5/48, CtlJ 1*0, CI21I80, C[3], 827/56, C 14] , *1 25/336 , 

C 15] 1*1/24, All], sl/5, A 12] I *2/5, 413], *1/3, Al4]|*4/5, A(5),«l, 

8 11,0] | b l /5 , B 12,0] ,80, 612,11(82/5, B 13 , 0] 1*7/36 , B13,1](*0, 

B 13,2] 1*5/36, Bl4,0],s0, B14, 11,80, Bl4,21,*4/5, B{4,3],*0, B15,0] 
,81/4, B 15, 1] ,*0, B 15,21 ,*»35/4, B 15 , 31 , *54/7 , B 15 , 4] , *25/ 1 4 , 



BUTCHER 2 2 


CC0]|a5/48| C[lJl«0» C[2J l«0r c 133 |«27/56* C 14] i«125/336» Ct51i = l/24 
» Atl]l*-l/5» A C2] |s2/5 f A C31 I « 1 /3 » A[4]184/5| At5]t«l» Btl,0]|3-1/ 

5 f Bl2,0]|a4/5 1 B[2,l] |s-2/5» B 13, 0] 1 87/36 > BC3,1]I*0» Bt3,2]*s5/ 

36 1 B (4# 0) I *0 > B[4,l]|s0j Bt4»2]1«4/5> Bt4,3]l«0» Bt5,0]J*l/4» 
Bl5,lll«0| B15,2 ] |b-35/4» B 15, 3] »b 54/7 f B t5*4] l«25/14» 

NEWTON-COTES, LAWSON 12 

‘ \ 

C 10] l«7/90f cm*«0f C 12] I *32/90 » C{3]l«12/90» Ct4]|«32/90» C15]I» 
7/901 ACt]lM/2f A t2] I * 1 /4 > A tS] t*l/2| A taS « «3/4 » At5]l*l> 

B Cl ,0] l»i/2f B [2,0] |83/I6f B t2,il I»l/I6f Bt3,01t*0l Bt3,l]|80j 
B [3,2] l a 1 /2 1 B (4, 0] t s 0 1 Bt4,l] |8-3/16f B 14 , 2] I a6/ 1 6 > B C4,3J l»9/16» 

B {5, 0} |sl/7y Bt5,l]»s«/7f Bt5,2]»s6/7| B [5,3! |8-12/7» BtS»4)l«8/7| 



NEWTON-COTES, LUTHER 13 


C [01 I *7/90 j C [ 1 ] |30| C (21 1 37/90 1 CC3}|»32/90* C[A]l«12/90| 

C [5] I *32/90 | A [ 1 ] | s 1 | A [2] t = 1 | A[3]|sl/4| A[4]|*i/2» A[5]|83/4f 
BC1#0] tal f B[2,0]i*l/2> B t2, 1] I s 1/2 f B [3,0] |8l4/64f B [3, 1] t *5/64) 

B [3,2] (8-3/64; B [4 , 0] l a- 12/96 1 B [« , i ] » 8- 1 2/96 1 B [4 , 2] I *8/96 j 
B [4, 3] (864/96) B15,0 ]|b 0» B [5, 1 ] ts-9/64 f B t5 , 2] I *5/64 | B[5,3]|s 
16/64 f B [5, 4] 1*36/641 

NEWTON-COTES, BUTCHER 2 

Cl0]»«7/90( C 1 1 ] I s0 f C [2] l*32/90| CI3]|al2/90f Ct4]|832/90f 
C [5] I »7/90 > All] l«i/8> A 12] I *1/4 f A[3]|*l/2( A[4]**3/4» A[5]|*l» 

B [1 ,0] 1*1/8 f B [2,0] i *0 1 B [2, 1) 1*1/4} B[3,0]i»l/2> B C 3 , 1 ] t 1 « 

B C3,2] jali B [4,0] |83/J6| B[4,l](s0; B[4,2]|*0» B [4 , 3] I *9/ 1 6 1 
B [5,0] fs-5/7 1 B [5, 1] l>4/7 f B [5,21 1*12/7 f B [5, 3] 18-12/7 l B[5,4]|*8/7| 



NYSTROM 


14 

C [0] 1*23/192, Ctl]|80| C (2) 1*125/192) C C31 « »0 9 C [4] •■- 81/1921 
C [5} 1*125/192) A(l]|*l/3f A 121 f «2/5 f AtS9l>l| A(4}t*2/3f A(5]t*4/5) 

B tl#0J 1*1/3) B [2,0] 1*4/25* B 12# 1 ] 1 *6/25f IB (3#0] 1*1/4) B (3* 1) )■> 12/4 1 
B [3,2] 1*15/4) B [4, 1 J 1*6/81 ) B Ca * l J 1 *90/81 ) B [4, 2} t*>50/81 1 B(< 4 , 3 ]|* 
8 /Blf B [5,0] 1*6/75) B [5, 1] 1*36/75 ) B 15,2] 1*10/75) B [5 , 3] 1*8/7 5 > 

B [5,4] | * 0 1 


SHANKSC6-6) 19 

C [0] 1*198/3696) C(U 1*0) C 12] 1*1225/3696.) C [3] I * 15 40/3696) 

C [4] 1*810/3696) C {5] (*-77/3696) C [6] i*0j C17]|*0) A [1] t *1/300) 
A C 2 ] 1*1/5) A [3] 1*3/5) A(4]|*14/15) A 15} 1 * i ) B 1 1 , 0] l *1 /300) 

B 12 , 0 ) l*-29/5) B 12, 11 1*30/5) B [3 , 0] 1*323/5 ) B [3# 1] t*~330/S) 

B (3*2] 1*10/5) 8 [4,0] 1**510104/810) B 14# 1] 1*521640/810) 

B (4,2] ts-12705/810) B (4, 3] I * 1 925/8 1 0 ) B [5 , 0] I *>4 1 7923/77 ) 

B (5, 1] 1*427350/77) B (5# 2] f *>10605/77 ) B (5, 3] t * 1 309/77 ) 

B [5,4] ) *>54/77 ) 



UT l 1 


C C 0 ] 1*124789/1000000) cm»«0» C 12] **574915/1000000) C 13] 1*1 1 1 13447/ 

10000001 C [4] ta-l0ei6403/1000000| C [5] l«3250/1000000| 

A (I] 13266000/1000000) A [2] 13398999/1000000) A (3) 13997475/1000000) 

A [4] »»1 ) A [5] 1*1/2) 

B Cl » 0] 13266000/1000000) B 12# 0] 1899749/1000000) 6 C2#l) 1*299250/ 
1000000) B C 3 » 0 1 1*997397/1000000) B C3, 11 I *-3740204/ 1 000000 ) fi[3,2!i* 
3740282/1000000) B C4#0] tal01855l/1000000 ) B(4,l) ) a-3826725/ 1 000000 ) 

B 14# 2] 133810768/1000000) B14#3] I *-2594/ 1 000000 ) 

BC5#0] t 3" 36199/1000000) B!5#11 1*956681/1000000) B[5#2] 13-519392/ 
1000000) B C5#3J *38376573/1000000) B [5# 4) I ••8277662/ 1 000000 ) 



1 


UT 2 

C [0] 1*78464/1000000, crn»« 0 | C [21 1*637651/1000000, C[3]|* 
93U08/1000000, C tai »»103537/1000000, C [5J i *-75076 1 / 1 000000 , 

A til I *203066/ 1000000, A 121 1*304600/1000000, A (31 1*620351/1000000, 

A 14] I • 1 1 A (5] t*l/2, 

Btl#0] 1*203066/1000000, B[2,0] 1*76150/1000000, 612,111*228450/ 
1000000, BI3,0]|*327625/1000000, B 13, 1 ] | *-10 1 6945/ 1000000 , 

B [3,2] 1*1309671/1000000, B t4 , 01 I *-300307/ 1000000, B [4, 1] t *2751 166/ 
1000000, B[4,2] (*-2664664/1000000, B [4, 3] I *1 2 13805/ 1 000000 , 

B *5,01 **293963/1000000, B [5, 11 t *-68779 1 / 1 000000 , B [5,21 1*891098/ 
1000000, 0 15,31 1*24998/1000000, B [5, 4] I *-22269/ 1 000000 , 


(T> 



UT 3 1 


C [0] 1*99466/1000000; CtlJl*0» C [21 I *3975 1 4/ 1 000000 ; C [31 | *242604/ 
10000001 C C4J t*-72772/1000000; C [5] 1*333186/1000000; 

A til 1*200933/1 000000 f 

A [2] 1*301400/1000000; A 131 l*606170/1000000f A(4]|*lf 
A [5] 1*918110/10000001 

Btl»01 »»200933/1000000| B (2* 0] I *75350/ 1 000000 f B [2, 11 1*226050/1000000 
I B [3*01 l*308207/1000000» B [3# 1] |»-934792/1000000; 

B C 3 # 2 ] 1*1232755/1000000) B [4*0] i*-232399/ 1 000000 ; 

Bt4*ll 1*2571643/1000000; B 14,21 I *-2609052/ 1 000000 f 
B [4*31 1*1269608/1000000; 

Bt5»0J 1*- 59480/1000000; B[5#ll 1*972686/1000000; B[5,2] la-663913/ 
1000000; B t5»3] 1*616247/1000000; B 15*41 1*52570/1000000; 


CT> 

ro 



SIXTH ORDER COEFFICIENTS 


FEHLBERG 6 5(6) 

C 1 0 ] «*7/1408f C 1 1 ] I *0 f C [2] 1*1 125/281 6| C[31*»9/32| C [41 | * 125/788 f 
C C 5 3 t *0 j C 161 1*5/66; C[7Jl*5/66f 

Atll I* 1 /6 

A [21 !■ « /15 » A [31 I* 2 /3 | A [41 |* 4 /5 , A [51 I* 1 ; 

A [61 |*0f A [71 |*1| 

BH,01 »• 1 /6 f B [2,01 is 4 /75 , B[2,il j« 16 /T5 f B[3,01i* 

5 /6 t B[3,1J i* -8 /3 f 013,21 I* 5 /2 | Bt4,0] i* *8 /5 ; 
B[4,1J |* 144 /25 f 014,21 |* -4 ; B[4,31 I* 18 /25 f 

B[5,0J 1* 361 /320 | B [5, 11 |* -18 /5 > B[5,21 I* 407 /128 » 
B[5,31. I* -11 /80 » B [5, 41 ** 55 /128 , 

B [6,01 l*-li/640» 

B [6, 1] 1*0; B [6,21 ;*1 1/256; B [6,31 l*-l 1/160; B [6,43 1*1 1/256; 

B [6, 51 1*0; B [7,01 l*93/640| B 17 , 11 I 18/5 » B [7 , 21 I *803/256 ; 

B [7,31 **-11/1601 B [7,41 »*99/256» S[7,51t*0; B[7,61;*i; 



FEHLBERG 6 6(7) 


c [0] 1=77/1440) cm 1 = 0 ) C[ 21 (= 0 ) C [3J 1 = 1771561 / 6289920 ) C[41(« 
32/105) C [5] 1=243/2560) C [6] l = i 6807/74880 ) C [7] 1*1 1/270; 

A tl] 1*2/33) At2Jisa/33» A[3](=2/ll) A[4](«l/2) A(5]{=2/3) A[6](=6/7) 
At7]»al| B[1,0] (=2/33) B[2,0]|=0) B(2,U(«4/33) B 13, 0J |»l/22> 
Bt3,i],=0, B[3, 211=3/22, B [4*0] 1*43/64? B[4*ll(=0) B [4,21 I.-U5/64, 

B [4*31 (=77/32) B [5*0] (=•2383/486) B[5*l]|=0f B [5*21 («1067/54) 

B [5*3] (=-26312/1701 ) B[5*4] (=2176/1701) B [6 * 0] | ■[ 0077/4802 ) 

B [6* 1] i =0 1 B [6*21 (=-5643/686) B [6* 3] 1 = 1 16259/16807) 

B[6*4] (■■6240/16807 ) B[6*5] 1 = 1053/2401) 8 t7*0] (■■733/176) 

B[7*ll(«0) B (7*21 ( = 141/8) B [7* 3] ) ■•335763/23296 ) B [7*4] (=216/77) 

B [7*51 (=•4617/2816) B [7*6] (=7203/9152) 



BUTCHER (6-7) 2 

C [0J l»ll/120f CtlJ*s0j C(21l»27/40| C(31l»27/40» C (41 |3»4/ 15 | 

C (51 I »■<♦/ 15 f C(6] |"11/120| Amssl/J| A (21 * *2/3 » A (31 t a 1 /3 J 
C (71 l«0> 

A(4)ial/2y A(5]l>l/2l A[61i«i» B 1 1 , 0 1 l s 1 / 3 > B(2,01t*0f B(2,llt*2/3» 

B [3,0J iel/12, B(3,llt*l/3> B (3,21 t*»l/12f B f 4 , 01 I *»1 / (6 f 

B[4, 11189/81 B [4,2J »*-3/16| B (4 , 31 t *-3/8 » B(5,03ia0* B(5,ll»*9/8» 

B (5,21 8.a-3/8| 8 [5 , 33 » *-3/4 » B (5, 41 I * 1 /2 » B (6 , 03 I *9/44 1 B(6,Ut« 

•9/ 1 1 1 B (6,23 I *63/44 } B [6, 31 8 * 18/1 1 f B(6,43s*0> B [6 , 5] 8 *«1 6/ 1 1 t 



BUTCHER(6-8) 2 


C [0] JS41/840, C til t*216/8«0r Cl2]|a0, C C3 J I *27/840 , C tfl] | *272/840 , 

C [5] 1827/840, C [6] 1*216/840 f C 171 t «4 1 /8a0 , AtlJ|«l/9» A(2) 181/6, 
At3Jl«l/3» A C4] ,»l/2, A [5] | s2/3 , Al6J|»5/6, At71l8l, B(l#0)|8t/9, 
B[2#01,8l/24, B I2» 11 t«3/24 , Bt3,0J|8l/6, B 13, 11 |s-3/6, BtJ, 21 184/6, 
B[4#01, 8 . 5 / 8 , BC4#l]|«27/8, B 14# 21 , b. 2 4/8 , B ( 4 * 3 ] 1 * 48/8 , 615 , 01,8 
221/9; B [5,1] ,8.901/9, B t5#21 1*867/9, B 15 # 3] ,8.102/9, B15#4J,*1/9j 
B C6#0J »s. 783/48, B (6# 11 ,*678/48, B (6 # 21 t 8.472/48 f B [6, 31 (••64/48 f 
B 16# 41 1800/48 , B [6#51 :b 3/48» B 17# 0] , *76 1 /82 , B 17 # 1 1 I s-2079/82 , 

BC7, 21 181002/82, B(7#3J 18834/02, B 17 , 41 1 8-454/82 , B [7 , 5 J , *-9/82 , 

B 17 # 6J 1872/82, 



SHANKS (6-6) 


19 

C [0] :*198/3696| cm 1*01 C (21 l»1225/3696f C (3] I * 1 5«0/ 5696 » 

C [41 1*810/3696) C 15] t*-77/3696> C(6]|*0f C(7]|*0| A(l]t*l/300y 


A(2]|*l/5y A 13] 1*3/5 f A [41 1*14/15 > A(5]l*l) B C 1 # 0 3 t*l/300f 
B (2,0] »*-29/5» 8(2,1] * *30/5 y B [3, 0] 1*323/5; B (3, 1) »»-330/5 y 
B (3,2] t*10/Sf B (4,0] l*-510104/8l0y B (4, 1J y*521640/810y 
B (4,2] y**»t2705/810y B (4, 3] |*1 925/8 10 y B (5 , 0] i *- 4 17923/77 y 
B (5, 1] l*427350/77| B (5,2] l*«10605/77y B (5,3] l*1309/77y 
B (5,41 »»-54/77y 



SHANKS (7"7) 19 


C 10] I ■ 1 4/300 j Ctl] H0| C £2] i*81/300f C £3] I * 1 1 0/300 f C£4]l*0> 

C [5] 1*81/3001 C [8] **14/300f Ct7]J*0> 

Atl] 1*1/192 » A [2] t * 1 /6 p A 1 3 J 1*1/2 p A[4]|*lp A£5]|*5/6p A£6]|*lp 
B ll»0] 1*1/192) B £2,0] »a«l5/6p B 12, 1] l*16/6p B [3*0] 1*4867/186} 

B 13, l] »*-5072/186p B £3,2] 1*298/186) B £4 , 0] I »• 1 9995/3 1 f 0l4,lJt* 
20896/31 p B £4,2] l*«1025/31 p B £4, 3] | * 1 55/3 1 p B (5, 0] I *-469805/5022 p 
B 1 5 # 1 ] 1*490960/5022; B 15 , 2] I *-22736/5022 p B 15,3] 1*5580/5022) 

B £5# 4] 1*186/5022 p B £6,0] I*9l4314/2604p B C 6 # 1 ] 1**955136/2604; 
B£6#2] |»47983/2604p B£6,3] t *«65l 0/2604 p B £6,4] I *<*558/2604 f 
B £6,5] l*2511/2604| 



SEVENTH ORDER COEFFICIENTS 


FEHLBERG 6 6(7) 

C 1-0] |all/064f Ctlj*a 0 ; Ct 2 ]*a 0 ; C [33 * e 1 77 1 56 1 /62S9920 ; C [4] * *32/ 1 05 ; 

C [5] *8243/2560; C [ 6 ] 1816807/74060 ; C C7I * 80 ; C C 8 ] »*1 1/270; Ct9Jj* 

1 1/270; 

All] ia2/33; A [2] 1*4/33; A[33*a2/li; A[4],*l/2, A[5]»a2/3; Al6]|86/7; 
A 17] 1*1; A 18] I 80 ; A( 9 ]| 8 l; 

011,01 182/33; B 12 , 01 * 80 ; B [2, 1] **4/33 ; B(3, 01**1/22; 

B [3, 1 ] * 80 ; 8 13,21 183/22; B [4,0] 1*43/64; Bl 4 ,l],a 0 ; B 14 , 2] ; a, 165/64 , 
Bt4,3] 1*77/32; 8 15,01 *8-2303/486; B15,1]*80, B 15 , 23 * a 1 067/54 ; 

B 15,3] *8-26312/1701 ; B 15,4] *82176/1701 ; B {6, 0J **10077/4802; 

B16, U*a0; B16, 2] ja-5643/686; B 16,3] *al 16259/16807; 

Bt6,4] *8-6240/16807; B 16,5] »«1053/2401; B 17, 0] * *-733/ 1 76 ; 

817,11*80; Bl7,23*ai4t/8; B 17, 33 » *-335763/23296 ; B 17 , 43 » «2 16/77 ; 
Bt7,5]*a-4617/2816; 817,61**7203/9152; 

Bt8,0] *815/352; BIS, 1]*80; B18,2]**0; 

B [8,3] J8-S44S/46592; B 18 , 4] * *1 8/77 ; B 18, 53 * a- j 2 i 5/5632 ; 

B [ 8 , 6 ] **1029/18304; B18,73*»0; B 19 , 0] * a- 1 833/352 ; B(9,l]*80; 

B [9, 2] 1*141/8; 8 19,33 *8-51237/3584; B 19 , 4] * = 1 8/7 » B [9 , 53 » «-729/5 12 , 

B [9,6] *81029/1408; B[9,73*a0; B[9, 83**1; 



FEHLBERG 6 7 (8 ) 

c [0] i = ai/aa0» cti]p*0p 01231=0, ct3ji*0p 0(431=0, c 1 53 1=34/105, 

C [6] I a9/35 , C (73 1=9/35, C (83 1=9/280, 0(9], =9/280, C [ 1 03 , =4 1 /840 , 

A 1 1 3 ! *2/27 , AC2]» = l/9f A(33l = i/6p A [43 I =5/ 1 2 , A C5J l = 1 /2 , A[6],a5/6, 
A [73 I * 1 /6 p A(8],=2/3p A[9]p = l/3p A[10]| = l, 8 C 1 * 0 3 I *2/27 p B(2,0]pa 
l/36p 8t2,l]|Bi/i2, B C3,03 l»l/2ap Bt3,lJ|e0p 8[3,2] l«l/8p 
B (4,0] pa5/12p B(4,i)l=0, B [4 , 23 * a-25/ 1 6 p B (4 , 33 I »25/ 1 6 p B[5,0]|a 
1/20P B [5, 13 l»0p B [5, 23 1 90 p B[S,3]«*l/4f 0(5,4], *1/5, B[6,0]|* 

■25/108 p B[6, 13l = 0p B [6,23 I ®0 , B [6, 3] I a 1 25/ 108 p B (6, «3 p a. 65/27 p 
816,53 »al25/5a, B [7, 13 l*3l/300p Bt7,lUa0, B[7,23l*0p 
B[7, 311=0, B[7,4J p*61/225p B 17,53 ia»2/9, 8 C7 , 63 t a 1 3/900 p 
B [8, 0J I »2 , 9 18, 13 l*0| B[8, 23 *80, B [8 , 33 I a-53/6 p B [8 , 4] | «704/45 , 

B 18,53 l«-107/9p B [8,6] ps67/90p B(8,73» = 3, B [9 , 0 J | *.9 1 / 1 08 p B[9,1 J,b 
0| B (9, 23 I »0 , 8 [9,3] pa23/108p 8 [9, 4] p a-976/ 1 35 p B [9,5] 1=311/54), 

B [9,6] pa-19/60, B [9,73 IS17/6P B [9,81 18-1/12, 8(10,0] 1 = 2383/4100, 
BC10, 131=0, 8(10,2] ! =0 , B ( 10, 3] pa-341/164, B [ 1 0 , 4 J , = 4496/ 1 025 , 

BUB, 5J pa. 301/82, B (1 0, fe] , .21 33/41 00 , B U 0, 7] » .45/82 p 0(10,0],, 
45/164, B [10,9] 1 = 18/41 p 



SHANK$(7-7) 19 

C (03 8*14/300) C(13l*0) C (23 1*81/300 j C 133 I * i 10/300 » C[43ia0f 
C (53 8*81/300) C [63 8*14/300) C(73 8*0| C 183 8 *0 ) Ct93l*0) C(i03)*0| 

A (13 8*1/192) A [23 8*1/6) A(33 8 *l/2j A[4]8*i) A[5]8*5/6) A[6]8*l) 

B (1,0] 8*1/192) B 12,03 (*-15/6) B [2, 13 8*16/6) B [3/0] 134867/186; 

B [3# 13 (*-5072/166) 6(3*23 1*298/186 f 8(4,03 )«-19995/31 | 8(4,138* 
20896/31) B (4,23 (*-1025/31 ) B [4,3] t>155/31 ) B [5, 03 1**469805/5022) 
B (5, 1] 1*490960/5022) B (5,23 (*-22736/5022) 8 [5,31 1*5580/5022) 

6 15,43 1*186/5022) 8 [6, 0] |*91 431 4/2604 ) B (6, 13 1**955136/2604) 

0 (6,23 1*47983/2604) B[6,33 I **65 1 0/2604 ) B (6, 43 8 *-558/2604 ) 

B [6,53 8*2511/2604) 



SHANKS (7*9) 19 


C t0] 1*1 10201/2140320, C [1 ] » 80 j C[2]|»0, C 1 3 ] ,*767936/2140320, 

C C4] 1*635040/2140320, C 15] *-59049/2140320, C [ 6 ] ,—59049/2140320, 

C t7] ,*635040/2140320, C 18] ml 10201/2140320| A(l)|.2/9, A[2),>l/3, 

A [3] ,*1/2, A [4] 1*1/6, A [5] ,*8/9, AI6]|*l/9f At7],*5/6, Atfl]|»lf 
B [ 1 » 0] | «2/9 » Bt2,0] mi/12, B [2, 1] t *3/12, Bt3,0],«l/8 , B(3,l],a0, 

B t3,2] 1*3/8 , B 14*0] , *23/21 6 , Bt4,l],*0, B [4 , 2] 1*21/216) 

B[4,3] |— 8/216, B[5,0] ,.-4136/729, Bt5,l],*0, B [5 , 2 ] » — 1 3584/729 , 

B [5/ 3] ,*5264/729, B t5,4] (*13104/729, B 16 , 0] i * 1 05 ] 3 i / J 5 1632 , 

B 1 ] 1*0, B16, 2] ,.302016/151632, B [ 6 , 3] , — 107744/151632, 

Bt6r4] ,—284256/151632, B 16,5] ,*1701/151632, B 17 , 0] , —775229/ 
13759201 817,13 1.0, B 17,2] 1—2770950/1375920, B 17,3] 1*1735136/ 
1375920, BI7, 4], .2547216/1375920, B C 7 , 5] 1*81891/1375920, Bt7,6],. 
328536/1375920, 0 18,0] ,.23569/251888, B[8,l],*0, B 16,2] 1 — 122304/ 
251888, B { 8 , 3 ] 1—20384/251888, B 18, 4] ,*695520/251888 , 

B [8,5] 1—99873/251888, B 18 , 6] , .*466560/25 1 888 , Bt8,7]|* 
241920/251888, 


— 1 
ro 



EIGHTH ORDER COEFFICIENTS 


FEHLBERG 6 7(8) 

C[0)l*0» C[1]|*0» C [2] * *0 1 Ct3]|*0| C[4]|*0» C 151 »*34/105j C[6]|* 
9/35» Ct7]|=9/35» C[8]l*9/280> CI9]|»9/280| Ctl01|*0| C 1 1 1 3 t *4 1 /840 f 
C [12] |*41/840> 

A [ 1 ] 1*2/27 f A [2] I * 1 /9 f A [3] I *1/6 | A[4]»*5/12j A[5]|*l/2f A[6]|*5/fe> 

A 1 7 1 1*1/6) A [8] I *2/3 f A(9]|*l/3f Atl03 »«1> A(ll]l*0f A C121 l a lf 

B 11,0] 1*2/27 | B [2,0] I* 

t/36» B [2, 1] «*1/12» B [3,0 ] *b 1/24> B[3,l]|*0, B[3,2]t*l/8> 

B [4,0] l*5/12| B[4,l)t*0| B [4,2] l*»25/16f B [4,3] 1*25/16* B[5,0]|b 
1/201 B [5, 1] 1*0} B [5,2] I *0 9 B [5,3] l*l/4f B(5,4]|*l/5> B[6,0]|b 
•25/108f B [6, 1] l*0| B [6,2] l*0f B [6,3] (>125/108) B[6,4) |B«65/27f 
B [6,5] |Bl25/5«| B [ 7, 1] 1*31/300 f B[7,l]|s0, Bt7,2]i=0» 

B [7, 3] |*0f B 17,4] |*61/225f 8 (7 » 5 ] I b* 2/9 | B (7 , 6] | * 1 3/900 f 
B [8,0] l*2f B [8, 1] |*0f B [8,2] I *0 f B [8 , 3] I *-53/6 f 8 [8 , 4] 1*704/45 , 



B[8,51 1 — 107/9, B [ 8 , 8 ] 1*67/90 ; B[8,71».3» B (9,0) (■.91/1081 B(9,1]|b 
0f BC 9 , 2 ]|* 0 ; 8 [9,3] |*23/108| B [9,4} |a.976/135f B (9, 5] I «3 1 1 / 54 , 

B t9 , 6] I a* 1 9/60 ( Bt9,7] 1*17/6; B [9,8] (.-1/12, B [ 1 0 , 0] 1 .2383/4 1 00 , 
BC10,l)|a0; B[10,2 J|b 0| B [ 1 0 , 3] I .*34 1 / 1 68 ( BU0,4] l*4498/1025| 

B[10, 5] ».-301/82> B[10, 8] 1*2133/4100; B [ 1 0 , 7] 1*45/62; Btl 0 , 8 ]|« 
«S/184; B [10,9] t ■ 1 8/41 ) 

B [1 1 ,0] f s3/205f B [ 1 1 , 1 ] | s 0 1 

B[ 11 , 2 ]|b 0 » B [ 1 1 , 3] 1 .0 1 B 111,4] 1*0; B [ 1 1 , 5] I *.6/41 ; 8 [ 1 1 , 6 ] | *.3/205 , 
B (1 1 , 7] »*»3/41 ; Bill, 8] 1*3/41; B [11,9] |*6/41 ; BC11,10]|*0; 

B [12,0] !* a 1777/4100( B[12,l]|*0| B[12,2]|*0; B 1 12, 31 t *«34 1 / 1 64 ; 

B [12,aj (*4496/1025; B112,5] (*.289/82; B [12,8] 1*2193/41001 
Btl2, 71 (*51/82; B [12,8] (*33/164; B [ 1 2, 9] » ■ 1 2 /a 1 ; B[12,10](*0; 

B CIS, 1 13 1*1; 



FEHLBERG* 8(9) 


C [0] ,*32256/1000000, CtU,*0, C [2] , *0 , C CS) I «0 f C[4],*0, C[5],*0, 

C[6],*0, C[7],*0, C 10] 1*259837/1000000) C [91 1*92847/1000000) 

C C103 I ■164523/1 000000| CC11J ,*176659/1000000, C [121 ,*23920 1/1 000000} 

C [131 * *3948/ 1000000, C 1141 ,*30726/1000000 , a til I "443689/ 1 000000 , 

A [21 J *665534/ 1000000, A [31 1*998301/1000000, A [41 I *315500/ 1000000, 

A [51 1*505441/1000000, A [61 1*171428/1000000, A [7] 1*828571/1000000 , 

A [8J 1*665439/1000000, A [91 1*248783/1000000, A [10] 1*109000/1000000, 

A [111 1*891000/1000000, A 112] 1*399500/1000000, A [13] 1*600500/1000000, 

AH4II*l, Btl,0]|*443689/1000000, B [2 , 0] | * 166383/ 1 000000 , 

B [2, 11 1*499150/1000000, 8 13,01 1*249575/1000000, B[3,l]f*0, 

B [3,21 1*748725/1000000, B [4, 0] t *2066 1 8/ 1 000000 , B[4,l],*0, 

B [4,2] 1*177078/1000000, B [4,31 ,*-68197/1000000, B [5,0] 1*109278/ 

1000000, B[5, 111*0, B (5, 2] 1*0, B [5,3] 1*4021/1000000, 8 15,41 ,* 

392141/1000000, B [6,01 ,*98899/1000000, B(6,l]|*0, B[6,2],*0, 

B [6, 3] 1*3513/1000000, B [6, 41 ,*124760/1000000 , B [6, 5] ,**55745/ 

1000000, B[7,0] ,*-368068/1000000, B[7,l],«0, B(7,2)t*0, B[7,31,»0, 

Cn 



B t7»a] **-2227389/1000000) B [7,5] **1 374290/ 1 000000 ) B [7, 6] 1*2049739/ 
1000000) B[8,0] **45467/1000000) B[8,l]t*0) B(8,2]**0) 8(8,3] |«0) 

B [8,4] 1*0) B [0,5] 1*325421/1000000) B [8 , 6] | *284766/ 1 000000 ) 

B [8,7] **9783/1000000) B [9, 0] * *60842/ 1 000000 ) B[9,l]|*0) 6{9,2]|*0) 

B [9, 3] **0) B [9, 4] 1*0) B [9,5] l*-21 184/1000000) B [9, 6] 1*193965/1000000 
» B [9,7] **-4274/1000000) B [9 , 8] ** 17434/ 1 000000 ) B (10,0] **54059/ 
1000000) B (10, 1] **0) B [10,2] **0) B (10,3] 1*0) B [ 1 0, 4] t *0 ) B(10,5]t*0) 

B [10,6] J*1 10298/ 1000000) B (10, 7] *«• 1256/ 1000000) 6 (10,8] 1*3679/ 
1000000) 8(10,9] l *-57780/ 1 000000 ) B [11,0] **127324/1000000* B(tl,l]** 
0) B (1 1,2] l«0) B (11, 3] *«0) BUI, 4] 1*0) B (1 1 ,5] 1*0) B(ll,6]|*0) 

B (11,7] **114488/1000000) B [11,8] **287730/1000000) 8 (11,9] **509453/ 
1000000) B (11, 10] **-147996/1000000) B [ 1 2 , 0] * *-3652/ 1 000000 ) 
B[12,l]t*0) B [12,2] *«0) B(12,3]**0) B(12,4]**0) B 112,5] 1*81629/ 



1000000, 8(12,6] (*-386077/1000000, B 112,71 I >30862/ 1 000000 , 

B (12,8] (*-58077/1000000, BC12,91 1*335986/1000000, 8(12,10] 1*410668/ 
1000000, B (12, 1 11 (*-l 1840/1000000, 8 (13,0] (■•1237535/1000000, 

B (13, 1] ,*0, B [ 1 3 , 21is0, B (13,31 (»0, 8(13, 4], *0, B (13,51 (*-24430768/ 
1000000, B (13,61 (*547795/1000000, B l 1 3, 71 ( *-4441 386/ 1 000000 , 

B (1 3, 8] 1*10013104/1000000, B (13,9) (s-14995773/1000000, B(13,10] ,* 
5894694/1000000, B(13,lll 1*1738037/1000000, 8(13,121 1*27512330/ 
1000000, 0(14,0] (*-352608/1000000, 8 (14,1] ,*0, 8(14,2)1*0, B114,3](* 
0| 8(14,3)1*0, 8(14,4] (*0, B (14,5] (*• 1 8396 1 / 1 000000 , B(14,6Ha 
• 655701/1000000, B 114,7] (*-390861/1000000, B ( 1 4 , 8] ( *267946/ 
1000000, 8 (14,9] (*-1038302/1000000, B ( 1 4, 1 0] (* 1667232/ 1000000 , 

B (14, 11) (*495519/1000000, B (14, 12J ,*1 139400/1000000, 8 (14, 13] (* 
51336/1000000, 
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C [0] 1*41/8401 CtlJ »B0> C [2] 1 «0 f C 13] t >27/640 f C [4] I *272/840 » 

C [5] l*27/840» C [8] »s2U/e«0> Ct7J»«0| CCS] |x21S/8a0r C C9J l*41/840f 
All] 1*4/27 f A C2] 1*2/9 f A(3]|*l/3j A[4]»sl/2 f A [5] t*2/3r A [6] (*1/6| 

A 17] (■! I ACS] I aS/6 i A[9]»*lf B t 1 , 0) | *4/27 | B t2 , 0] » *1 / 1 8 ( 

B 12# 1 J 1*3/18} 8 [3,0] 1*1/12 f B[3,l]|s0 f B [3, 2] 8 *3/ 1 2 } B [4 , 0] ( *1 /8 f 
B[4,l]|*0( 8ta#2]»»0f B[4,3]»*3/8( B [5,0] |*13/54| B[5,l]t*0f Bt5,2] 
l*-27/54( B[5,3J 1*42/54} B [5, 4J 1*8/54} B [6, 0] i *389/4320 } B[6,l]t*0} 
B [6,2] 1**54/4320} B 18,3] 1*966/43201 B {6# 4] 1 8*824/4320 1 
B t6#5] l*243/4320f B [7,0] t«*23l/20} B[7,1J|*0} B [7,2] t*Bl/20f 
B [7# 3] 1 8*1 164/20 } B [7,4] 1*656/20} B 17,5] 1**122/20? B(7, 6)1*40} 

B 18,0] 18-127/2881 B C8 , 1)1*0} B 18, 2] I *18/288 ) B C 6 , 3 1 18*678/288; 

B [8,4] |845b/288f B 18,5] l*-9/28B| B [0, 6] 1*576/288 f B 18, 7] 1*4/288 | 

B [9,0] 1*1481/820} B [9, 1) 180) B (9, 2] (■•81/6201 B [9,3] t*7104/820» 

B [9,4] l*»3376/820} B 19,5] 1*72/820 } B 19, 6] » **5040/820 } B[9,73l»*60/ 
820 ( B [9,8] l*720/820f 
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C [0] I &41/840 * C(1]**0, C (21 * *0 , C (31 **0* C(4]**0* C (5] **216/840* 

C C6] **272/840* C I7J **27/840* C 181 ! *27/840 * C (9] * *36/840 * C1101I* 
180/840* C till **41/840* C 1 121 1*0* A Cll l*l/9| At2H«l/6| At3]|*l/4, 

A 14J 1*1/10* A 15] **1/6* A (6] 1*1/2* At71**2/3* A[8]**l/3| A(91**5/6, 

A 1 1 01 * *5/6 * A [111**1* B (1 , 0] **1/9* B 12*01 1*1/24* B (2# 1) 1*3/24* 

B 13, 0] **1/16* B(3, 1]**0> 6(3,21**3/16* B (4, 0] * *29/500 * 8t4,l]**0* 

B [4,21 **33/500* B 14,31 ***12/500* B (5#0J **33/972* B[5, 11**0* 

BIS, 21**0* Bt5,3J **a/972* B 15,41 1*125/972* B t6,0] *«• 21/36, 

B 16, 11 **0 * BC6, 21**0, B 16,31 **76/36* B t6 , 4] * • 1 25/36 » 8(6,51** 
-162/36* 817,0]* *.30/243* B(7,i),«0* B(7, 21**0, B (7, 3] **-32/243, 

B (7,41 **125/2U3* Bt7, 51**0* 8(7,61**99/243* 0(8,03**1175/324* 

B(8, 1],*0, B (6 , 2} * *0 * B (8,31 **-3456/324, B (8 , 4] * *-6250/324 * 

B ( 8 , 5 ] **8424/324* B (8 , 6 ] » *242/324 * B (8 , 71 * *-27/324 * B (9 , 0] * *293/324 * 
B(9, 11**0* B(9, 2]**0* B(9, 3] *«- 052/324* B (9,4] **-1375/324, 



B t9,Sl l«1836/32ay B [9 , 6] t*»118/324* B [9 , 71 1*162/324* B C9 , 83 I «1 f 
B 110,0] 1*1303/1620* B[10,lll«0> B(10,2]**0* B [ 1 0, 3] I **4260/ 1 620 * 

B 110,4] **-6875/1620* B [ 1 0, 51 * *9990/ 1 620 » 8 1 1 0, 6] I * 1030/ 1 620 * 
Btl0,7]»*0> B [ 1 0 # 8 ] * =0 * B 110,91 1 = 162/16201 B ! 1 1 , 01 | *-8595/4428 » 
Btll, 11**0* Bill, 21 |b 0| B 111,31 l = 30720/««28» Btll, 4] 1 = 46750/4428* 
B tli,5] >*-66096/4428* Bill, 61 **378/4426* B C 1 1 , 7} (••729/4428) 
Btll,8J >*-1944/4428* Btll, 9] **-1296/4428* Bill, 10] 1*3240/4428* 

CURTIS 3 

Ct01l«l/20* Ctll to0* C (21 1=0* C [3] t «0 * C [4] * «0 * C (5] 1*13/160 f 
C [6] *«36/180f C [71 1 = 64/180* Ct81l«l/5| C [9] I * 1 3/ 1 80 * Ctl0]|«l/20| 
C c 1 1 3 t «0 1 C 112] 1=0* A 111 1*183855/10000001 A [2] 1*183855/1000000 * 

A [31 >*275775/1000000* A [41 1*689439/1000000* A [5] 1*827326/1000000* 
A [6] 1*172673/1000000* A [71 1*1/2* A [8] 1*827326/1000000* 



A [ 9 ] |al72673/1000000| A C10] * = 1» B 1 1 , 0] I a 1 83850/ i 000000 ) 

B I2»03 »*92128/1000000| B [2, 11 1*91922/1000000) B [3,0] (■68943/1000000) 
B[3»l] 1*0) 8 [3,2] (■206831/1000000) Bt4»0] ( *689439/ 1 000000 ) 8(4,11 i* 
0) B [4,2] (*-2585396/1000000) B 14 , 3] ) *2585396/ 1 000000 ) B (5 , 0] ( *82732/ 
1000000) B [5 , 1 ] ( *0 ) B (5 , 2] ( *0 ) B E5, 3] (=413663/1000000) B(5,4](* 
330930/1000000) B 16, 0] (*97115/1000000) B(6,l])*0) Bt6,2](*0) 

B 16, 3] (*97308/1000000) B [6,41 (*-44005/1000000) B (6, 5] ( *22254/ 
1000000) B [7,0] (*-623200/1000000) Bt7,l](*0) BI7,2](*0) 

B [7,3] (*-259948/ 1000000) B [7,4] 1*249582/1000000) B 17,5] (*-113452/ 
1000000) B [7,6] 1*686138/1000000) B [8 , 0] ( *308 1 86/ 1 000000 ) 8 [@,1]|*0) 

B [8,2] 1*0) B [8,3] (*519705/1000000) B t8, 4] I *-76 1 906/ 1 000000 ) 

B[8,5] 1*419861/1000000) 8 [8,6] (*-629640/ 1000000) B[8,7] 1*971118/ 
1000000) B[9,0] (*139493/1000000) B[9,l](*0) B[9,2](*0) B[9,3]|* 
117240/1000000) B [9, 4] (*-249427/1 000000) B[9,5] (*-26481/1000000) 

B [9,6] I *-11 8354/ 1000000) B [9,7J 1*182542/1000000) B [9,8] (*127657/ 
1000000) B [10,0] (*-499040/1000000) B[10,l](*0) B[10,2](*0f Bll0,3](* 
• 1386394/1000000) B [10,4] (*1331109/1000000) B [ 1 0 , 5] 1 *-674026/ 
1000000) B110, 6] 1*1119619/1000000) B 110,7] (*-592592/1000000) 

B [10,8] 1*506298/1000000) B (10,9) 1*1195027/1000000) 


00 
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